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(1) Introduction. 

The .purpose of this report is to introduce the reader to the prin- 
ciples of operation and capabilities of birefringent filter systems. 
Birefringent filters have a number of significant advantages over more 
conventional spectral filtering techniques based upon gratings, Michelson 
interferometers, Fabry-Perot interferometric filters, or all deposited 
interference filters. Their primary advantage is the wide field of view 
compared to the other interference filters. The small angle sensitivity 
yields throughput gains of between 50 and 300 for imaging or collecting 
light from extended sources. This means a factor of 50 to 300 in number 
of resolution elements or total light collected in the desired spectral 
band. Besides the throughput advantage, birefringent filters have the 
capability of being tuned over 1.25 decades in frequency, can have 
tailored spectral profiles , and can have a set of selectable bandwidths . 
The devices are virtually insensitive to the direction of gravity. 
Furthermore, the filters can easily operate from 0°C to 50°C without any 
active temperature control. 

Birefringent filters have long been used by solar physicists to 
study spectral characteristics of the solar atmosphere, -mostly in the 
Hydrogen alpha line. Their use by scientists in other disciplines, 
however, has been very limited for a number of reasons. Commercial 
filters have been normally available only for a single wavelength., 
although two commercial tunable filters have been made. The filters 
are expensive, with single— frequency , 25 ram clear aperture filters 
costing about $50,000. Perhaps most importantly, the operation of the 
filters has not been understood by potential users and the basic sources 
of information about the filters have required knowledge of wave propaga- 
tion of polarized light in isotropic media at a discouragingly high level. 
In addition, the basic work on birefringent filters has been published 
in foreign language journals for which only private translations are 
available . 

We feel that the major limitation on the use of birefringent filters 
has been the lack of a single reference where an interested scientist can 
obtain an overall picture of the principles of operation and capabilities 
of these systems. This report is written to provide this overall view. 
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/ '. Normally, birefringent filters are explained using the techniques for 
analyzing propagation of polarized light, namely, the Jones Calculus. 

In this report the Jones techniques are avoided because of their 
unfamiliarity to the general reader and, perhaps more importantly, 
because of the lack of insight they provide into the behavior of bire- 
fringent plates and polarizers. While the Jones Calculus does an 
excellent job of determining the action of polarizing systems, it 
often does not make it physically clear why the optical system 
operates as it does. Fortunately, birefringent filters can also be 
explained by means of pulse response techniques, which only require 
some bookkeeping and a limited familiarity with the ideas of Fourier 
analysis. The result of this approach is a better physical under- 
standing, but less mathematical analysis capability. 

Chapter 2 is devoted to explaining the three fundamental types of 
birefringent filters : the Lyot (perfect internal polarizers) , the Sole 
(no internal polarizers), and the filter with partial internal polarizers. 
The most straightforward of the filters, the Lyot, is used to introduce 
the basic concepts and analysis methods for birefringent systems. The 
filter with no internal polarizers is discussed next. Here it is 
expected that the reader will have to make some effort to follow the 
somewhat tedious, but reasonably coherent explanation. Chapter 2 ends 
with a discussion of partial polarizing filters which, while mathemati- 
cally more complex than Sole filters, represent a small logical step 
from the Sole. 

Chapter 3 is concerned with how the bandwidth or full width at 
half m axi m u m (FWHM) , free spectral range (FSR) , finesse (gF) , and profile 
shape are determined for birefringent filters. The main point of the 
chapter is to give the reader the tools to determine the basic design 
parameters of a filter system. Explored in some detail is the subject 
of profile shape which is a parameter which is not usually tinder the 
control of the filter designer, whereas for birefringent filters, 
virtually any profile shape can be achieved. 

Chapter 4 discusses the subject of filter tuning. The methods for 
tuning the three basic filter types are established and some design 
implications of various tuning techniques are discussed. It is assumed 
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in the tuning section that achromatic tuning components are available. 
Specifically, it is the ready availability of easily fabricated achro- 
matic wave plates xrfiich makes the methods discussed in this section 
possible. 

The birefringent widefield element is introduced in Chapter 5. 

While a thorough analysis requires a large amount of tedious algebra, 
a simplified summary is presented as briefly as possible. The details 
of the field of view effects are deferred to an appendix, while only the 
important results and design considerations are presented in Chapter 5. 

Chapter 6 introduces the subject of bandpass control. Here it is 
demonstrated that it is possible to construct Lyot and partial polarizing 
filters which can have a preselected set of bandwidths. Sole filters 
are not discussed because freedom to adjust all plates means that the 
bandpass of a Sole filter can be continuously adjusted. 

The engineering considerations involved in reducing the theory to a 
workable filter system are discussed in Chapter 7. While this chapter 
does not describe all the considerations well enough, to enable a novice 
to build birefringent filters , it does indicate potential difficulties 
and their solutions. Since birefringent filters are sophisticated 
optical instruments they require some care 1 in their fabrication. Chapter 

r« 

» - ( , 

7 will indicate, however, that7their tolerance requirements can be much. 

less severe than- other, 'comparably powerful, '"optical systems. 

Chapter 8 provides a concise summary of the capabilities of bire- 
fringent filters.- It highlights their operating principles and advantages 
over other filter systems. In addition, some of the current work in 
extending the range of applicability of these filters is outlined. * 

Finally, two appendices are included which discuss in more detail 
some topics which are felt to be important. Appendix 1 is a detailed 
analysis of the field of view properties of birefringent filters. 

Although the results of this analysis are summarized in Chapter 5, this 
is such an important aspect of these filters that its inclusion is bene- 
ficial. The engineering details, specifically for space qualifiable 
filters, are discussed in Appendix 2. There, a more thorough presenta- 
tion is given than was included in Chapter 7 . 
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Throughout this report, spectral transmission characteristics are 
always plotted for a full free spectral range of the shortest element in 
the filter. The plots are always for transmission with-respect to phase 
retardation in the shortest element, which corresponds directly to the 
frequency of the incident radiation. 

This report, then, should yield understanding of the physical 
properties and control techniques for birefringent filters. With this 
knowledge the reader should be able to specify but not design a filter. 

For those desiring -more detailed knowledge of the details of filter 
operation and analysis, we .have prepared a companion volume to this 
report that represents a reasonably complete set of scientific papers 
on birefringent filters. It includes translations of most of the 
original Lyot and Sole work, and reprints of many of the important 
recent papers. Space only permitted the reproduction of papers speci- 
fically devoted to filter techniques. For background material on 
polarized light and the Jones Calculus techniques, the excellent summary 
by Shurcliff and the reprint collection edited by Swindell are recommended. 
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(2-A) Lyot Filters (Perfect Internal Polarizers) 

Birefringent filters were introduced by Lyot in 1933, and also 
developed independently by Ohman in 1938. These filters have been 
described in the literature, particularly in the papers by Lyot (1944) 
and Evans (1949) which are reproduced in Volume XI of this report. The 
basic theory is repeated here because of its simplicity and as an intro- 
duction to the pulse response technique of birefringent filter analysis. 

A Lyot filter is made up of a sequence of elements , each of which 
consists of a birefringent plate between perfect polarizers. As will 
be true of all birefringent plates in this discussion, the plate is 
assumed to be a uniaxial crystal cut with its optic axis in the plane 
of its faces. For a Lyot element the plate is placed with its optic 
axis oriented at 45° to the orientation of the parallel entrance and 
exit polarizers. When viewed through, a spectrograph the transmission 
of such an element appears as a series of alternate light and dark bands, 
a channel spectrum. Specifically, the transmission, t, of the element is: 

t(X) = cos (=-“> (!) 

with A = (n o - n & ) , C2) 

•wfrase d is the thickness of the plate, a is the difference between the 
ordinary, n Q , and extraordinary, n g , indices of refraction, and X is the 
wavelength. A thick plate produces a channel spectrum of many thin bands 
closely spaced while a thin plate produces a channel spectrum of fewer 
thick bands more widely separated. 

A Lyot filter, in its simplest form is made of a cascade of n such 
elements, each twice the thickness of its predecessor. The important 
feature of a cascade of Lyot elements is that the perfect polarizers 
isolate the elements so that the total transmission function is simply 
the product of the transmission functions of all the elements. Thus, if 
the elements have thicknesses 

d fc = 2 k d o k = 0, n - 1 , (3) 
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the total transmission function becomes 


n-1 Tr2 k d A 

t(X) = II cos^C —-) . (4) 

k=o A 


For maximum transmission at an operating wavelength, the total 
retardation, d Q , of the thinnest plate should be an integral number of 
wavelengths. The thickness, d Q is chosen to satisfy 


d 

o 



(5) 


for some integer, m. All thicker plates will automatically satisfy the 
same criterion. 

Shown in Figure 1 is the transmission function of each of the elements 
of a four— element Lyot filter, together with the cumulative transmission 
from each subfilter of the cascade* These are shown in decreasing order 
of thickness of the elements. It is clear from the figure that the 
action of each successively thinner element is to double the separation 
between transmission maxima while changing the shape of each maximum only 
slightly . 

The separation between transmission maxima is the free spectral range 
(FSR) of the filter. As is evident from Figure 1, the FSR of the overall 
filter coincides with the FSR of the thinnest element. Substituting (5) 
into (1) we see that at the operating wavelength, the transmission 
function is 


rAd 

T 


r (A ) = cos^ ( — r — — ) = cos^Cmir) = 1 . 


The next transmission maximum will occur at a wavelength 


C6) 


\ - X o + X 


fsr 


5 


^Remember that all transmission functions are with respect to phase retarda- 
tion in the thinnest element, which corresponds directly to frequency. 
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such that 


irAd 

o 

A +A , 
o rsr 


= (m-l)ir. 


In usual practice 


X » x c 

o fsr 


and (.7) may be approximated by: 


irAd X 

(1 ^~) = (m-l)TT , 

A o o 


C7) 


which implies that the free spectral range, X becomes 

fsr’ 


fsr Ad 


C8) 


If it is assumed that the profile of the transmission maximum at A q 
is the same as that given by the thickest element then the filter's bandpass 
may be calculated. If the half width at half maximum of the transmission 
profile of the thickest element is A , then (1) , (3) , and (5) imply: 


and that 


2 

cos 


7r2 n "^Ad \ 
o hw / 


1 

2 


n-1 


Ad 


X +X, 
o hw 



(9) 


( 10 ) 


Again assuming that X >> A, 

o hw 


A. 

hw 



(ID 
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The filter bandwidth, or more precisely, the full width at half maximum 
(FWHM) is 


FWHM = 


2 M i 


( 12 ) 


In the subsequent analysis this estimate will be refined to account for 
the entire filter instead of just the thickest element, but (12) is a 
good approximation. Using (8) and (12), the finesse, 3F, of a Lyot filter 
may be approximated as 


_ *n 

FWHM 


(13) 


In practice the best Fabry-Perot interferometric filters with apertures 
greater than a few millimeters have a finesse of 20 to 30 . This is 
comparable to a 4 to 5 element Lyot cascade. 

As mentioned in the introduction, a technique which has proven to 
be especially useful in optical filter analysis is impulse response or 
time domain analysis . This is a technique used regularly by electrical 
engineers and geophysicists, but only recently by optical designers. 

Impulse response analysis is complementary to the more common spectral 
analysis, although the two are mathematically equivalent. Either the 
time domain or the frequency domain response is sufficient to characterize 
a system since the two are conjugates under the Fourier transform operation. 
The freedom to investigate a filter's behavior in either domain is often 
quite useful. In particular, much of the experience of electrical 
engineers in time domain analysis is directly applicable to birefringent 
filter analysis and synthesis. Surprisingly, impulse response analysis 
of birefringent filters is actually much simpler than impulse response of 
more general electrical networks. 

In Figure 2 the impulse response of a single Lyot element is diagrammed. 
A unit amplitude impulse (indicated by an arrow) emerging from the entrance 
polarizer has two components, one along and one orthogonal to the optic 
axis of the plate. These each have amplitude 1/-/2. Light propagates in 
these two polarizations with velocities c/n Q and c/n^ respectively. As a 
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Figure 2. Pulse propagation through a simple Lyot element. 
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result the two component impulses incident on the crystal emerge at 
different times with a separation. At, given by 



(14) 


When these two impulses encounter the exit polarizer, only their components 
in the pass polarization orientation are transmitted. The transmitted 
components again have relative amplitudes The response of the entire 

element to the single incident unit impulse is therefore, a pair of pulses 
of amplitude 1/2 separated in time by At as given by (1.4). A thick 
element produces pulse pairs with long separations while a thin element 
produces closely spaced pairs. 

When two Lyot elements, whose length ratio is 2:1, are cascaded, 
the output from the first element serves as input to the second. One 
input pulse creates four output pulses. More generally, the n element 
cascade discussed above will have 2 n output pulses for each input pulse. 

In Figure 3 the pulse response for a four-element Lyot filter has been 
illustrated. There, as is often done for convenience, the time axis has 
been redefined such that the output pulses axe symmetric about the time 
origin. Half the output thus leads the input while the other half lags. 
This is a mathematical convenience and has no physical significance. 

It is clear from Figure 3 that the factor of 2 length reduction between 
successive crystals is the property which creates equally spaced, identical 
output pulses at each stage of the cascade. 

Some properties of the impulse response are evident from Figure 3. 

The spacing of the pulses is 



05) 


and is determined by the thinnest element. The total time duration, T, of the 
response is 


T = (2 n -liA't . 

o 


06 ) 
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Figure 3. Pulse tree of a four-element Lyot filter. The pulse output of each element 
of the cascade is shown. 
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When, n is moderately large so that 
2 n » 1 

Then T is approximately 

T “ 2t ' (17) 

where t is the time delay of the longest element. 

In both the spectral analysis illustrated in Figure 1, and the pulse 
response shown in Figure 3, it was assumed that the longest element was 
first, followed by successively shorter elements. It should be evident 
that the spectral response, being the product of single element spectral 
responses, is unchanged by a rearrangement of elements. Equally true, 
however, though perhaps not obvious, is the invariance of the pulse 
response to rearrangements . The detailed structure of the pulse response 
tree of Figure 3 changes, but the final pulse output will remain unchanged. 
It is a useful exercise to verify this for at least one rearrangement. 

As mentioned previously, the pulse response and the spectral response 
are intimately related to each other through the Fourier transform opera- 
tion. Since we have already derived the two responses for Lyot filters, 
we can demonstrate this connection. Some mathematical notation and 
definitions will be needed and are summarized below. 

An impulse at time t k will be represented by the Dirac delta function, 
6 (_t - t^) , which has the properties : 


6 (t-t k ) - 0 for t y t^ 


/ f (.t) 6 Ct-t k ) dt = f ( t k ) . (18) 

—00 
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The Fourier transform, F(v), of a function, f(t), is defined as: 


oo 

FCv) - /ict) 

— CO 

with the inverse transform being: 


(19 a) 


03 ' 

f U) = 2 “ / *60 e“ i2irvt dv . 

—CO 


(19b) 


Any specific properties of the Fourier transform will be introduced as 
needed. 

Using the above definitions, the equivalence of the pulse and spectral 
responses of a single Lyot element may be demonstrated easily. Assume that 
the pulse response of Figure 2 is represented by 

fCO = ^K t+f% +|d(t (20) 

This represents a pair of pulses separated by At as given by (14) The 
Fourier transform of (20) is 


F(v) 


CO 

1 f,, . At. “i2irvt, . 

- J 6 (t + y-) e dt + 

—CO 

1 -iirvAt , 1 iirvAt 

tt e + r e 


1 f At. “i2 , n'vt J 

2 J 6 (t “ 2 ~ ) e dt 

—co 


- COSirvAt . 


( 21 ) 
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F(v) as given by (21) is the amplitude response of the element. The trans- 
mission is given by the squared magnitude of the amplitude response which 
becomes 

2 tfd(n-0 

t(A) = cos (22) 

c d 

when we substitute, v = t-> and. At = — (n -n ) . This is the same trans- 

a c o e 

mission function as given previously, Cl), for a single Lyot element. 

For a two-element Lyot filter the transform of the pulse response 
is: 


FCv) = \ I 


6(t + | At) + SCt + | At). 


+ «(£ - | At). + 6Ct - | At). e i27fVt dt 


1 „ . 1 
= — cos 3irvAt + — cos irvAt 


(23) 


where At is again given by (14) , and is the time delay of the shorter 
element. Using the trigonometric identity 


. . _ „ ,A+B. ,A— Bv 

cos A + cos B = 2 cos (—^V cos (— — ) 


(24) 


and the relation for v the transmission becomes 


/ ^ 2 2ircAt 2 ircAt 

t(v) = cos — 7 cos 


X ' 


(.25) 
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The equations (21) and (23) suggest that for an n element Lyot 
cascade there will be a sum of 2 n impulses in the pulse response, a sum 
of 2 n cosine -terms in the amplitude spectral response, and a product 
of n cosine squared factors in the transmission function. Though this 
may be proven by repeated use of identity (24) , here only the result 
will be presented. 


,n-l 


l n 1 k=l 


cos 


(2k~l)-?rcAt 

X 


n-1 

II 

k=o 


cos 


„k 

2 -rrcAt 


(26) 


It is not true in general that any sum of cosines can be written as a 
single product, but the relationship between these sum and product forms 
will be used in further arguments-. The specific result, (26), for the 
Lyot filter is a consequence, again, of the 2:1 length ratio between 
successive elements. 

Using the pulse response —spectral response relationship, it is easy 
to find a closed form expression for the product of cosines, - (26). . In 
this case it is convenient to redefine the time origin to coincide with 
the earliest output pulse, resulting in a transformed response: 

00 Z -1 

' F(v) 6 (t-kAt) e i2irvt dt 

2 k=o 

2 n -l 

_ 1 i2irkvAt 

- ~*r Z e c 

2 k-o 

Since this is simply the sum of a geometric progression it has the closed 
form evaluation 
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F(v) 


■1 , i2ir2 n vAt 

1 1 - e 

„n i2irvAt 

2 1 - e 


1 sin2 n 7rvAt 
2 ^ sinirvAt 


iir2 n vAt 
e 

iirvAt 

e 


As the transmission is simply the squared magnitude. 


(28) 


TOO 



. 2 ,2 ircAt x 
sm ( — 7 ) 


2/ircAt( 
sm (— — ) 


(29) 


Using expression (29) it is possible to re-derive some of our earlier 
results, but with more accuracy. In particular, ,it is evident that unit 
intensity transmission -maxima occur whenever is an integer. This is 
the same result as (.6), Now, however, it may be seen that the half 
maximum transmission occurs wh.eneyer 


„ cAtA, 

. 2 ,„n hw N 
sm (2 u — r ) 

1 „ 1 *1 

2 ~ , n cAtX, 

4 . 2 , h^ 

sxn ( r — ) 

*1 


(31) 


When 2 >>1 it is possible to simplify (31) further, yielding 


cAtA, „n cAtA, 

/oil. hw s 2 hw 

sm C2 ir — •) = 


a2 

A o 


/2 H 


(32) 


This is a transcendental equation of the form 


sxn q = q 

/2 


( 33 ) 
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whose solution is 

q = 1*39 (34) 

This implies that the full width at half max i minn is 

*2 

FWHM = .44 j (35) 

2 cAt 

which is a more accurate result than the earlier estimate, (12) , based on 
the length of the longest element alone. 

The free spectral range is , as given previously in (S) , 


fsr cAt 


Consequently the finesse is now given as 


& = 1.13 2 


n 


(36) 


The subject of Lyot filters and their modifications will be returned 
to in later chapters. 
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(2-B) Sole Filters (No Internal Polarizers) 

. The Lyot filter discussed previously was invented in 1933 and was 
developed into a practical solar instrument during the period 1933 to 
1950. In 1953 the next major advance in birefringent filter design 
came when Sole realized that a filter could be made without internal 
polarizers. The explanation of the basic operation of the Sole filter 
first appeared in English in 1958 by Evans. By 1964 Harris and his 
associates had developed a general procedure for the synthesis of 
arbitrary spectral response profiles with Sole-type filters. In this 
report all such filters without internal polarizers will continue to 
be called Sole or Sole-type filters even though their development owes 
much to Harris and others whose contributions will be described sub- 
sequently. As will be seen, the principle of operation of a Sole filter 
is far less intuitive than that of a Lyot filter. This is reflected in 
the literature on the subject prior to 1967, none of which provided a 
reasonably clear explanation of the mechanism. Fredga and Hogbom (1967) 
introduced the concept of pulse paths and path amplitudes which provides 
a simple qualitative explanation of the Sole filter and a reasonably close 
quantitative approximation for most Sole designs. Volume II of this 
report contains reprints of translations of the original Sole papers 
from the Czechoslovakian Journal of Physics, in addition to a Sole paper 
which appeared in JOSA, the Evans paper mentioned earlier, the Harris 
synthesis procedure, and the 'Predga and Hogbom paper. 

The original Sole filters, introduced in 1953, were made as a cascade 
of identical birefringent plates, at specified orientations, between 
entrance and exit polarizers. In this discussion, the angular orienta- 
tions will be denoted as: 


<f> o - orientation of entrance polarizer 
(usually taken to be zero) 

<j>k - orientation of the fast axis of the 
k-th plate 


* 


P 


- orientation of exit polarizer 


Sole presented two different types of filters which', while equivalent, 
are not identical. These are known as the fan and folded filters and 
are specified as: 
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Fan filter — 4* Q = 0 

4 k = (2k-l)a l<k<n (1) 

<j> = 0 

P 

Folded filter - «j> = 0 

<i> k = (-l) k a l<k<n (2) 

♦p " 90 ° 

where a = 45°/n (3) 

for both configurations. 

In Figure 1 the transmission profiles for several fan type filters 
are shown along with the corresponding spectral responses of equal length 
Lyot filters. The figure demonstrates some interesting properties of fan 
filters. The first is that they are not dramatically different from the 
Lyot filters. This is a surprising feature to many who are accustomed 
to the Lyot design with internal polarizers. The second property is that 
for very few elements the sidelobe level is lower than the corresponding 
Lyot profile, but for moderate and large values of n the sidelobe level 
is higher. This behavior was noticed early in the study of these filters, 
and the design was criticized for its relatively poor sidelobe behavior. 

As will be shown later, the reason for this behavior is not too difficult 
to understand and the mechanism to improve the performance follows 
easily. 

Transmission profiles for the folded filters are not presented 
since they have a remarkable equivalence to the fan filter profiles. A 
folded filter has the exact same transmission profile as a fan filter 
except that the spectrum is shifted by one half of the free spectral range. 
For this reason, fan type filters will be the principal subject of this 
discussion. 

Evans (1958) presented a general formula for the transmission profile 
of the fan (and thus the folded) Sole filters. This transmission is: 
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SI3SL08BJ 



PISTES 
RET VHG 

1 I. 00 45.00 

2 2.00 45.00 
POLRRIZERS 

RHO RNG 

1 0.00 0.00 


90 135 180 225 270 3tS 360 


9.0 SIOELOae* 0.S2E-01 



PLATES 
RET RNG 

1 1.00 45.00 

2 2.00 45.00 

3 4.00 45.00 
POLARIZERS 

RHO RNG 

1 0.00 0.00 

2 0.00 0.00 


90 135 180 225 270 315 380 


18.0 SIDEL06E* 0.48E-01 



PLATES 
RET RNG 

1 l. 00 45.00 

2 2.00 45.00 

3 4.00 45.00 

4 0.00 45.00 

paimuzERS 

RHO RNG 

1 0.00 0.00 

2 0.00 O.QQ 

3 0.00 0.00 


45 90 135 ISO 225 270 315 380 


3, 5, and 15 element fan Sole 
.ement Lyot filters (d,e,f) . 
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T(X) » f sin n X 


l 


sxn x 


cos x. tan 


•]' 


where 


cos X = cos a cos 


TTCAt 


(4) 


and At is the time delay of a single birefringent plate as in the Lyot 
filter analysis. Transmission maxima occur, as stated earlier, at 
wavelengths for which cAt = kA. Formula (4) was derived using a power- 
ful technique known as the Jones Calculus (h. C. Jones, (1941)), but, 
while elegant, the derivation gave us little insight into the mechanism of 
the Sole filter. 

The lack of an intuitive explanation was typical of the literature 
on Sole filters. While the Jones Calculus was a powerful analytic tool 
it did not indicate the enormous potential of the Sole filters. 

Sole (1959) showed that by a slight modification of the angles 
given by (1)— (3), the sidelobe level could be reduced by a factor of 
5 or more. This provided an answer to the earlier criticism of the filters, 
but only touched on the real capability, namely that by a suitable modifi- 
cation it is possible to achieve almost any spectral profile. 

As examples of the capabilities of Sole filters, three responses 
have been selected. The first response, shown in Figure 2(a) is the 
transmission profile of a seven element Sole filter which is equivalent 
to a three-element Lyot filter. Both have identical, uniform pulse 
responses and, consequently, identical transmission profiles. Figure 2(b) 
is the transmission profile of a filter whose pulse response has a linear 
taper, or triangular shape. With increasing time delay the impulse 
response rises linearly to a central maximum and then decays linearly 
back down again. This results in a further sidelobe reduction over the 
uniform pulse response. If the pulse response amplitudes are made pro- 
portional to binomial coefficients, (^)-> a transmission profile with no 
sidelobes at all results. This profile is shown in 'Figure. 2 tc)-. In order 
to facilitate comparisons, the orientations and selected profile charac- 
teristics are reproduced here: 
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Basic Sole 

Lyot 

Triangle 

Binomial 

♦l 

6.4° 

8.4° 

3.5° 

0.6° 

*2 

19.3° 

18.9° 

11.8° 

5.0° 


. 32.1° 

31.4° 

25.7° 

19.5° 

*4 

45.0° 

45.0° 

45.0° 

45.0° 

*5 

57.9° 

58.6° 

64.3° 

70.5° 

*6 

70.7° 

71.1° 

78.2° 

85.0° 

*7 

83.6° 

81.6° 

86.5° 

89.4° 

FWHM/FSR 

.113 

.112 

.144 

.199 

Sidelobe 

.089 

.052 

.003 

0.0 


As the figure and the table show there is no limit to the sidelobe reduc- 
tion which is possible with suitable angular modifications from the basic 
Sole design • The implication is that with the appropriate procedure for 
determining the plate orientations, any transmission profile is available. 
In the following^ two techniques for determining plate angles corresponding 
to a given response will be described. 

An elementary description of the Sole filter operation was provided 
by Fredga and Hogbom (1967). Here, for the first time, the concept of 
pulse paths and path amplitudes was introduced. Harris et al (1964) had 
previously used pulse response techniques in their synthesis procedure, 
but not in an intuitively appealing development. Unfortunately, the 
simplicity of the pulse response analysis of Lyot filters doesn't quite 
transfer to Sole filters. The difficulties arise because the plates 
are not isolated from each other by polarizers as they are in the Lyot 
configuration. This means that at the interface between two plates the 
impulse transmission from one plate to the other must be considered in 
two orthogonal polarization states instead of just one. Mathematically 
this means that the impulse response must be considered as a vector 
quantity instead of just a scalar. 

In this section each of the two polarization components will be 
considered separately. Any pulse incident on a plate will be decomposed 
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into two pulses , one along the fast axis and one along the slow axis of 
the plate. This is the natural coordinate system for the plate and is 
the key to an understanding of the entire system. The logical consequence 
of this is that a pulse incident on the first plate produces two output 
pulses. These in turn result in four pulses output from the second 
plate, eight from the third, and finally 2 n from the last plate and 
from the exit polarizer. They, of course, occur at only n+1 distinct 
time delays so there is considerable overlap. It is much easier, however, 
to calculate the amplitude of each of the 2 n pulses individually than to 
try to calculate the pulse response directly. This is the path amplitude 
approach, since each of the 2 n output pulses can be imagined to be the 
result of a specific path through the system. Each pulse will have 
followed a path through some fast axes and some slow axes. The paths 
through a two-plate filter can be diagrammed as : 


f h — 1 — I 
sill 


f t- 
s \ 


f I H fill 

s l — 1 S 1 — M 


V f 2 


V s 2 


s l :f 2 


s l ls 2 


(5) 


The advantage of this approach is the ease of calculating the pulse ampli- 
tude associated with each path. This is illustrated by a too -plate 
example. 

In Figure 3, a general Sole configuration with two plates between 
entrance and exit polarizers is illustrated. The respective orientations 
are 4> q , ^ . In the analysis of such a system the absolute orientations 

are unnecessary since the filter must be invariant to overall rotations. 

The quantities of interest are the angle differences which will be denoted; 


0 k = " *k-i 

e = ^ - a . 

p p n 


( 6 ) 
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Figure 3. Pulse propagation through a general two-element Sole filter. 
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With this notation, the response of the filter to an incident impulse 
may be considered. Such an impulse will generate two components incident 
on the first birefringent plate; one along the fast axis direction and 
one along the slow axis direction. These components will have relative 
amplitudes : 


AC^) = cos 0 1 
A(s^) = sin 0^ 

respectively. Each of the two pulse outputs of the first plate may be 
similarly decomposed into components along the fast and slow axes of the 
second plate and the resulting amplitudes calculated. These become: 


A ^ f l :f 2^ = cos S 2 

A(Si : f 2 ) =-sin ©2 

ACfi :s 2 ) = sin © 2 
ACs 1 :s 2 >- = cos ©2 


A(f 1 ) = cos © 2 
ACs^) =-sin ©2 
A(f x ). = sin © 2 
A(s^X = cos' ©2 


cos ©^ 

sln0 i 
cos ©^ 

sin ©^ . 


( 8 ) 


Finally the exit polarizer at orientation, 
components whose amplitudes are: 


ij> , will transmit impulse 


A(f^:f 2 :p) = cos ©^ A(f^:f 2 ) 
AG s i:f 2 :p) = cos Q p A(s 1 :f 2 ) 
ACf-^:s 2 :p) =-sin 0^ ACf^:s 2 ) 
ACs^:s 2 :p) =-sin 0^ A(s^:s 2 ) 


= cos 

© 

P 

COS 

9 2 

COS 

e i 

--cos 

6 P 

sin 

0 2 

sin 

°i 

=-sin 

© 

P 

sin 

0 2 

cos 

e i 

=-sin 

0 

P 

cos 

0 2 

sin 

V 


The notation used in (7) - (9), the colon notation, indicates the 
local orientation of a scalar impulse component as it goes through each 
plate. This is the pulse path concept introduced by Fredga and Kogbom 
in which, the output impulse response is decomposed into an enumeration 
of all possible paths through the plate system. With, each path, is 
associated an amplitude, as in (9)-> and a time of arxiyal. The time of 
arrival depends on the number of slow axes trayersed by the. path, since 
each, slow axis traversal results in a relative delay. At, with, respect 
to a fast axis traversal. In keeping with, the convention that the 
impulse response of a system is centered at the time origin, the output 
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time for any path will be determined by the simple rule: 


Rule 1: Each, fast axis traversal contributes a 

time delay of -At/2 

Each, slow axis traversal contributes a 
time delay of +A't/2 , 

In addition, the notation, D(f Q ), will be used to mean 5(.t-t o ), an 
impulse at time t Q . The total impulse response for a general two-plate 
system becomes, therefore: 


A~(t). = cos 0. cos 0„ cos 0 D(rAt). 
z ± z p 

-[sin 0, sin 0„ cos 0 

X JL p 


+ cos 0. sin 0„ sin 0 ] D(.0X 
1 J P 

i 

P 


ao) 


1 ^2 p- 

- sin 01 cos 0 . sin 0 D(iAt) 
-L Z P 


Since the spectral response depends only on the relative amplitudes of the 
pulse response, it is convenient to normalize (10) such that the earliest 
output pulse has unit amplitude. The response is then: 


A 2 (.t) = DC=-At) 

- .[tan © 1 tan 0 2 


+ tan 0 O tan 0 J D(0). 

2 p 

- tan 0 n tan 0 DCfAtX. 
1 P 


Cll) 


It should be apparent from (11) and the preceding development that 
the pulse response may be derived by some simple rules for evaluating 
path amplitudes simply from the plate to plate transitions of each path. 
The form of these rules will be taken to be the normalized form, to 
produce the amplitudes (11) rather than (10).. The specific rule is: 


Rule 2: The first plate contributes a factor 1 to 

the fast path and a factor tan 0^ to the slow 
path 

Each subsequent transition contributes a factor 
according to: 
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f 

I 

f 

s 


i-l :s i 


s . i :s 4 
i-l 1 


1 

tan 0. 
i 

-tan 0 # 

l 

1 


The exit polarizer contributes a factor according to : 


s :p -tan 0 . 

n p 

finally, the normalization constant is: 

n 

cos 0 II cos 0. . 

P i i 

Using these rules the total pulse response may be computed by a brute 
force enumeration procedure# Though this is tedious, the pulse response 
may often be approximated without a total enumeration, and as a theoreti- 
cal tool, the path enumeration procedure has been very valuable. Many 
general properties of Sole filters are proven most easily By arguments 
concerning individual path amplitudes. As the enumeration is conceptual 
rather than computational in those cases, the efficiency is not an issue. 

For an n element Sole filter there are 2 n possible paths but only 
n+1 possible time delays . In order to understand the pulse response from 
the path enumeration it is necessary to divide the set of paths into 
subsets of equal time delay# 

The path which only traverses fast axes appears at the relative 

time, - 7* At. Similarly the all-slow path, emerges at time, + ~ At. These 

z 2 z 
are, of course, unique. The output pulse at time, At, will be the 

resultant of the paths which include exactly one slow axis# There are 

n of these since there are n plates. In general, the output pulse at 

relative time - ~ ~ - At k = 0,...,n will Be the resultant of (£) 

different paths. 

Now, if Rule 2 is considered, it is evident that the f:f and s:s 
transitions do not attenuate the relative amplitude, hut that f :s and s:f 
transitions contribute a tan0 attenuation. For the usual case in which ©<45° 
this is indeed an attenuation. This means, that for each time delay, those 
paths with the fewest axis type changes (f:s or s:f) will contribute the 
largest amplitudes to the resultant output pulse. 
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To first order, therefore, the paths which dominate the pulse response, 
in increasing order of time .delay, are giyen as: 


A(f 1 :f 2 : . . 

• :f :p) = 1 
n r 

t 

11 a . 

= - 2 At 


A (f^ i f £ • ■ • 

.:f , :s :p) 

n— 1 n r 





= -tan© tan0 

t 

- -2=2 it 



n p 


2 


A(f ^:f 2? . • 

. :f , :s . :s :p) 
n-k n-k+1 n r 





= -tan© . , - tan© 

t 

n-2k A _ 
— At 



n-k+1 p 


2 


A(f ^5 « « » 

.:s :p) = tan0^.tan0 
n 2 p 

t 

_ n-2 . 

- — it 


A(si : s 2 : .. 

. :s :p) = -tan0-tan0 
n Ip 

t 

XI . 

= 2 At. 

(12) 


A first order approximation to a filter design for a given pulse 
response would result from the n+1 uncoupled simultaneous equations (12) . 

Take, for example, the construction of a seven-plate filter whose pulse response 
is equivalent to a three-element Lyot filter, that is, eight equal 
amplitude pulses. Equations (12) imply that the ©^ are all equal and 
that 


(13) 

is maximized, so 


(14) 

and 

0 =0.-90 

P i 

= -78°. 75 - 0-5) 

This is the filter whose entrance and exit polarizers are parallel at 
0° orientation and the seven plates are uniformly spaced at angles 

$ k = 11°. 25 k k = 1, ... ,7 . (16) 


1 = - tan 0. taxi 0 # 
i P * 


This is satisfied, and, since the normalization constant 
.is the transmission, when 

o =-22 

i n+1 


= 11 . 25 
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Since the design criteria, equations (12)., were only approocimate.it is 
not surprising that this filter only approximates a uniform amplitude 
pulse response. In particular the pulse response of this filter is: 


Time (At) 

75311357 
.2 2 2 ~ 2 2 2 2 2 

Amplitude 

.17 .13 .11 .09 .09 .11 .13 .17 


The spectral response of this filter is shown in Figure 4. As 
indicated there, the FWHM is 10.1% of the free spectral range, and the 
highest sidelobe is 13% of the peak transmission. This is contrasted 
with the spectral response of the three element Lyot filter shown in' 
Figure 2(a). There the FWHM is 11.2% of the free spectral range but 
the sidelobe level is only 5.2% of the peak transmission. 

A closer look at the impulse response of this first order approxi- 
mation is quite revealing. The pulse amplitudes taper downward from the 
outside inward. This means that the central pulses are being attenuated 
much more than the outside pulses even though their dominant paths are 
all equal. The reason for this is that equation (12) has accounted for 
paths of the form 


f 

s 

f :s 

but not paths of the form 


f 

s 


f 

s 


"1 I — i f 

J — — I J 

f :s s :f f :s 


(17) 


(18) 
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Figure 4. Transmission versus frequency for 
plate distribution. 
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a seven-plate Sole filter with a uniform 
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since these paths are of third degree in tan©^ instead of first degree* 
There are, however, many such paths for the interior pulses (remember 
that there are paths for the k-th pulse) and they are all negative 

since they have one s:f transition* In particular there are 

* 

(n - k + 1) (k - 1) (19) 

paths of the form (18) for the* k~th pulse* Since this is a quadratic 
function of k it gives the first order approximation its characteristic 
form* This pulse response is unfortunate since it is just the opposite 
of the response which minimizes sidelobes* Responses which minimize 
sidelobes peak in the center and taper downward towards the ends* In 
order to make a good filter, therefore, the angles must be modified to 
compensate for this effect* The Sole fan filter does this by making the 

« 

first angle difference, 0^, only half the others (while adjusting all the 
angles to maintain symmetry of course) . 

It is possible to incorporate the type (18) path amplitudes into a 
set of equations similar to (12) and solve for the higher order approxima- 
tion to a given pulse response. This isn’t too difficult, but it is a 
little too cumbersome to include here. Another, related, technique is 
to start with the first order approximation and to iterate to achieve a 
desired pulse response. This is possible since, for small angles, 0., 
the angular corrections, A0^ , which need be made at each step are approxi- 
mately linear functions of the pulse response error terms. Thus, at least 
for filters which have no extreme deviations from the basic Sole response, 
the. Sole technique of making small angle modifications will converge to 
the desired filter. This is the result which Fredga and Hogbom present as 
a consequence of the path amplitude arguments they introduced: sidelobe 

reduction to any desired level may be achieved by small angle modifications 
from the basic Sole design, or, more generally, any n-HL term Fourier series 
approximation to an- arbitrary transmission profile may be achieved by 
iteration from the first or second order filter configurations. 

In 1964 Harris, Amraann, and Chang went further, but in a way which 
is less intuitive. They used some reasonably sophisticated mathematical 
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techniques to relate the impulse response to the roots of n-th degree 
polynomial equations . One result of this is a rigorous mathematical proof 
that any of the above mentioned transmission profiles can be realized by 
a Sole- type filter. Another result, however, is that for most trans- 
mission profiles the filter is not unique. There are essentially differ- 
ent plate angle configurations which give tlie same transmission profile. 
This result is non- intuitive and not widely realized, but has some 
practical implications for actual filter construction. 

Consider a filter constructed to give a specified 1 transmission 
profile, [CQOI . This filter will have an impulse response whose 
pulse amplitudes are: 


C i i = 0, . . . ,n . (20) 

If the exit polarizer were rotated 90° the filter would now have a 
transmission profile which is the complement of the original, 

|r»CA) I 2 1- |c(A)| 2 (21) 

and an impulse response whose pulse amplitudes are denoted 

D i i - 0, . . . ,n . (22) 

.... . o 

The essence of the argument is that the square magnitude, | C (x) J » or, 

| , does not uniquely specify the impulse responses, or D^, 
respectively. The filter is, however, specified by the impulse responses. 
Consequently any of the compatible and Ih sequences will generate a 
different filter with the specified transmission profiles. In Figure 5(a), 
the transmission profile for a 9-element filter with a uniform pulse 
response is presented along with two different orientation configurations 
which achieve that profile. Notice that one of ■the angle sequences, <fi^, 
shows the monotonic progression typical of a modified Sole configuration. 
The other, however, has larger angle differences and a much less regular 
pattern. Both of these configurations imply the same transmission profile, 
|c(X)! 2 , rejection profile, jD(X)j^, and transmission impulse responce, C^. 
They differ, however, in their rejection impulse responses, D^. The 
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Figure 5* Transmission versus frequency of a nine-element Sole filter with a uniform pulse 
response (a) and its seven-element approximation (b) shown on linear (a,b) and 
logarithmic (c,d). scales. 
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impulse responses for these two filters are giyen as : 



♦i 

6.5 

14.4 

23.7 

34.1 

45.0 

55.9 

66.3 

75.6 

83.5 


I 

c. 

X 

.10 

.10 

.10 

.10 

.10 

.10 

.10 

.10 

.10 

.10 


D i 

-.01 

-.03 

-.04 

-.06 

-.09 

-.11 

-.14 

-.18 

-.22 

.88 


>1 

15.5 

47.9 

77.4 

77.1 

45. 

12.9 

12.6 

42.1 

74.5 


II 

G. 

X 

.10 

.10 

.10 

.10 

.10 

.10 

.10 

.10 

.10 

.10 


D. 

X 

-.03 

-.10 

-.19 

-.26 

-.18 

-.09 

-.43 

-.40 

-.53 

.36 


The existence of these multiple solutions has some practical impor- 
tance. Some of the solutions may turn out to be more, tolerant of alignment 
errors than others, and this question is a current research topic. Another 
important question raised is whether any of the solutions have consecutive 
plates oriented at the same or nearly the same angle. This has the practical 
significance that a double length plate could be substituted for the two 
single length plates. This would save on fabrication complexity. In 
Figure 5(b), such a seven-element approximation to one of the nine-plate 
filters of Figure 5(a) is shown. Notice that two pairs of internal plates 
have been merged into double length plates at their average orientations. 

The performance is nearly identical. 

The above example is not an unusual occurrence. In Figure 6(a) two 
filters which achieve the same, triangular, impulse response are presented. 
Again, there are plate pairs which are very nearly parallel and may be 
combined as shown in Figure 6(b). Again, the performance 'is nearly identical, 
with the seven-plate filter actually doing a little better. 

Another important aspect of the Harris synthesis procedure is its 
value as a diagnostic tool. Since it determines the plate angles for a 
given transmission profile it may be used to adjust a filter. If the 
filter's transmission profile is measured, the synthesis algorithm may be 
used to determine the angle at which each plate is currently set. It is 
a simple matter, then, to readjust the angles to their desired settings. 
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Figure 6. Transmission versus frequency of a nine~element Sole with a linear taper pulse 
response (a) and Its seven— element approximation (b) shown on linear (a,b) and 
logarithmic scales (c,d). 
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(2-C) Partial Polarizer filter 

The effects of using imperfect, or partial, polarizers in birefringent 
filters were discovered in 1953, the same year as the introduction of the 
Sole filter, Giovannelli and Jefferies discovered in an analysis of Lyot 
filters that in some cases the secondary maxima were suppressed if 'the 
internal polarizers were imperfect (the ratio of transmissions in the 
rejection and pass directions was non-zero) . In 1965 Beckers and 
Dunn verified the partial polarizer result via a computer program for 
Lyot filters. However, experimental attempts to verify the effect 
failed and the result was generally considered a curiosity. In 1974 
Title gave a simple theoretical explanation of why a single partial 
polarizer worked in a Lyot filter, and later developed several filter 
designs based upon multiple partial polarizers. A number of filters 
were constructed using partial polarizers and the theoretical predic- 
tions were verified. Previous experiments were shown to have failed 
because of neglect of the birefringence of the sheet partial polarizer. 

A partial polarizer differs from a perfect polarizer in that it 
not only has transmission in the pass direction, say the X direction, 
but also transmits in the orthogonal, Y, direction. Mathematically, 
the partial polarizer has two transmission amplitude coefficients, p 

X 

and p y . Usually P x is near. unity and p y is much less. The key to 
understanding the action of a partial polarizer is to consider it 
equivalent to the parallel combination of two optical paths. The first 
path will have a perfect polarizer with transmission amplitudes 


i 



p y = o . a) 

The second path will have a neutral density filter with, a transmission 
amplitude 


I 


P “ 
n 


P 


y 


C2) 


along all directions. 
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A general description of the effect of partial polarizers follows 
directly from the pulse path enumeration technique deyeloped in the 
previous section. Each partial polarizer will have an easily describable 
effect on each of the pulse paths. Before developing the general 
analytic technique, however, it is useful to consider an important 
special case; the two-element Lyot filter with its intermediate 
polarizer replaced by a partial polarizer. 

In Figure 1 the two-element Lyot filter with partial polarizer is 
shown in its conceptually equivalent form. It is depicted as a pure 
two-element Lyot in parallel with a pair of plates separated by a 
neutral density filter. 

In a normal Lyot filter the plates may be at either plus or minus 
45° with respect to the polarizers . Since the perfect polarizers isolate 
each element from its neighbors the relative orientation of successive 
plates does not affect the filter. When the polarizers' are partial, 
however, the angles between the plates become important and must be 
considered in the analysis. This means that for the two-element example 
two cases must be distinguished; the plates are either parallel or 
perpendicular. 

In the path of Figure 1 which contains the neutral density filter 
the time delay will be the sum of the time delays if the plates are 
parallel but will be the difference if the plates are perpendicular. 

The impulse response is, therefore, 

A® = ^ [ DC _ 3|t) + DC _ |tj + BC |t } + D( 34t } j 

+ ft [d(:- At) + D ] (3) 

where the sum and difference in the second term corresponds to parallel 
and perpendicular respectively. The transmission of the filter is then. 
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Figure 1. Optical schematic of a two-element partial polarizing filter and its conceptually 
equivalent form of a parallel two-element Lyot and two elements with a neutral 
density filter between them. 
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x ,2-rrcAt, /irCAt. 
- P ) cos C— — )- cos 


+ p cos (■ 

y 


(2±l)ireAt 

X 


>]• 


(4) 


Although (3) shows the distinction between the perfect polarizer and 
neutral density filter components of the pulse response, it is not a 
form which is very convenient for analysis. If, however, (3) is 
rewritten. 


A„Ct) = [ (1+2S) D(- ~At) + D(- ^At) + D(|At) 

+ (1+2S) DcfAt)] (5) 


Aj_(t) = C^) [D(- |At) + (1+2S) DC- ^At) 

+ C1+2S) D(^At) + D(-|At)] (6) 


with 




P -P 

x y 


it is obvious how the transmission will be related to the polarizers. 
Since p^, the rejection direction transmission amplitude, is always less 
than p^, the pass direction amplitude, S is always positive. The 
parallel plate pulse response, (5), has a taper which is high at the 
ends and low in the center. As in the previous section, this taper 
produces large sidelobes. The perpendicular plate response, (6), has a 
taper which produces lower sidelobes. Specifically it is higher in the 
center than at the ends. 

In Figure 2 the transmission profiles for parallel and perpendicular 
configurations are shown for several values of S. As may be seen, the 
parallel plate transmission sidelobes. increase with S as expected, while 
the perpendicular plate sidelobes decrease. In fact, when S equals one, 
the sidelobes disappear entirely. The pulse distribution is proportional 
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Figure 2. Transmission versus frequency for a partial polarizing module with 
crossed (a,b,c,d) and parallel plates (e,f,g,h) for S=0.0, 0.111, 
0.5, and 1.0 (p=0.0, .1, .333, and .5). 
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to 1:3;3;1, a binomial distribution, and binomially distributed pulse 
responses bare no sidelobes. For this minimum, the polarizer trans- 
mission amplitudes are in the ratio 



and the transmission intensity ratio is 



(7) 


( 8 ) 


which is quite a poor polarizer. 

A temptation, after seeing how much sidelobe control is available 
in a two-element Lyot configuration, is to replace more polarizers with 
partial polarizers in a multi— element Lyot. It is possible to use 
partial polarizers in a three-element Lyot filter and to achieve sidelobe 
reduction, but for more elements the use of partial polarizers raises the 
secondary maxima. There are designs, however, in which some of the polar- 
izers are partial and some are perfect, P y = 0, which are beneficial. 
Alternate partial polarizer configurations which are discussed in Title 
(1974), are very attractive. 

The discussion above has been concerned with the special case of 
Lyot configurations. More generally it is possible to consider Sole 
configurations with the possibility of partial polarizers between some 
of the plates. Since it is known that the Sole configuration -already 
has enough degrees of freedom to allow the synthesis of arbitrary trans- 
mission profiles, the addition of the polarizers would only increase the 
multiplicity of solutions for any given profile. This feature, however, 
might make it possible to find filter configurations which minimize the 



i ■ r* 
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In order to analyze the general plate/partial polarizer . configura- 

i 

tions, the pulse path, technique of the previous- section has been useful. 
The partial polarizers introduce some more notation, but actually very 
little more complexity into the procedure. The notation introduced in 
the previous section is: 

<J> o - orientation of entrance polarizer 
- orientation of k-th plate l<k<n 
$ - orientation of exit polarizer 

To this must be added: 

~ orientation of the k-th partial polarizer 
(following k-th plate) l<k<n-l 

Q P x~ P 

i, ' C ) for the k-th partial polarizer 

x y 

Notice that a perfect polarizer is represented by Q = 1 while no polarizing 
is represented by Q = 0. For general partial polarizing, 0<Q<1. With 
this notation, the rules for constructing pulse path amplitudes become: 

A(i v p k : W = c ° sC *k+i-V + \ cosC Wr 2 V 
ACf k :p k :s k+i ) ’ sin< WV + Q k sln «k^kfr 2 V 
A( WW + Q k sln « k + * fct . 1 -2* k ) 

ACs k :p k :s k+i ) = cos »k+rV - Q k cosC *k + i'kfr 2 *k ) e> 
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and as before: 

ACp:^)- = cosC4> 1 -^ 0 ) 

A(p:s 1 ) = sin(^> 1 -$ o ) 

ACf^ip) = cosC^ p -$ n ) 

A(s n :p). =*-sin($ -$ n ) . (10) 

Equations (9). and (10) represent the amplitude factors corresponding to 
each transition. To find the amplitude of any particular path the 
corresponding factors must be multiplied together. 

While it is possible to derive equations (9) and (10) from first 
principles, similar to the derivation of the Sole filter rules, they are 
also directly accessible from the Sole results themselves . A plate- 
partial polarizer— plate combination is essentially the same as a three- 
plate Sole configuration. The middle plate, however, has zero retardation 
introduces the additional attenuations , p and p , for paths along its 

x y 

respective axes. Since no time separation is introduced by the middle 
plate, the two paths which would normally he present are easily combined, 
resulting in the four path, amplitudes, (9), instead of the eight which 
might otherwise be expected. 

Equation s ^9) simplify considerably for some special cases. For 
example, in the perpendicular Lyot configuration discussed earlier the 
Q k angles alternate between +45 and —45 . In that case, if the partial 
polarizers are all oriented at 0°, the terms of equations (9) simplify 
such that each factor has only one term: 

A(f :p:f) = Q k 
A(f:p:s) - 1 
A(s:p:f) =-l 

A(s:p:s). =-Q k . ( 11 ) 

The particular case for which — 0 for all polarizers in the system 
is, actually, a general configuration. In this case, the transmission is 
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a maximum since, when the plate retardations are multiples of the wavelength 
the plates axe effectively absent. The transmission is only a function of 
the polarizers for those wavelengths, and is -maximized when all polarizers 
are parallel. 

In general partial polarizer designs are important, both theoretically 
and practically. The extra degrees of freedom may allow for some easily 
implementahle realizations of useful transmission profiles. Additionally, 
relaxation of the perfect polarizer requirement, p =0, usually re sul ts 
in a significant increase of p^. This produces dramatic increases in the 
overall transmission efficiency of the filters.- 

At the present time a general synthesis procedure is being developed 
for partial polarizing filters. The extra degrees of' freedom add some 
complexity over the Sole filter synthesis since they must be constrained 
either by restrictive assumptions or by optimization criteria. Both 
approaches are being used, and a systematic investigation into partial 
polarizing configurations will soon Be possible. 
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(3) Transmission Profile Management 

In the previous chapter a theoretical foundation for the three filter 
types has been presented. While some of the design criteria were mentioned, 
their implications and control were largely ignored. Some of these criteria, 
specifically bandwidth, free spectral range, secondary maxima, and overall 
transmission, will be discussed in this chapter. Other performance objec- 
tives such as the field of view and' tunability will be deferred to subsequent 
chapters . 

Of all the performance criteria, the easiest to control is the free 
spectral range (FSB.) . In the frequency domain the response is a periodic 
function, which is directly related to the impulse train response in the 
time domain. The periodicity in the frequency domain is, therefore, the 
reciprocal of the pulse position spacing in the time domain. Usually, the 
pulse position spacing will be the minimum distance between pulses in the 
impulse response. In general, however, the pulse position spacing will be 
the greatest common divisor of the interpulse spacings. If this greatest 
common divisor is At, then the corresponding free spectral range, FSR, will 
be 


FSR 


cAt * 


a> 


Since there is a direct linear relationship between pulse spacings and 
crystal lengths, the time delay, At, will correspond to the greatest common 
divisor, of all the crystal lengths. Specifically 


At 



d 

o 


c 


(21 


For a standard filter configuration all the crystal lengths are multiples 
of the shortest length.’ Most frequently then, the free spectral range is 
simply determined by the length, d^, of the shortest crystal according to 
(1) and (2). 

In order to increase the free spectral range of a filter it is necessary 
to make the shortest crystal even shorter. For a Lyot or partial polarizing- 
configuration this usually means adding one more short element to the filter. 
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For a Sole, no internal polarizer, filter with identical plates, each 
crystal must be made shorter and the entire design reconfigured. Composite 
designs, in which a thin Lyot element is used as a blocking filter to 
increase the free spectral range of a Sole filter, are also feasible. 

The design criterion which is usually of greatest concern is the 
bandwidth, or full width at half maximum (FWHM) . Although bandwidth is 
intimately linked to secondary maxima, or sidelobe, level, it may be approxi- 
mated separately . While the free spectral range is determined by the minimum 
interpulse spacing, the bandwidth depends mainly on the total time duration 
of the impulse response. 

The minimum bandwidth response for a given total time duration, T, 
results from a single Lyot element of length 


d 


T 


cT 



For this element the transmission is 


(3) 


T(lj = cos 2 (^) 


and the full width at half maximum is 


FWHM = 


2cT 


(4) 


05) 


where \ is the operating wavelength. 

As the crystal of length, d , is divided into shorter crystals to make 
a Lyot, Sole, or partial polarizing filter, the FWHM will increase somewhat 
while the secondary maxima recede. If the total length, d^, i S used to make 
an n element Lyot filter, the bandwidth becomes approximately: 

™Lyot " X - 17 ^ ■ <« 

Similarly, if the total length is sliced into n equal length elements and 
used for' a Sole fan or folded filter, the bandwidth becomes: 
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FWHM e _ » 1.59 

Sole 


2cT 


(7) 


In general, except for extreme impulse response apodizations , the FWHM will 
be inversely related to the total crystal length and will be of the above 
form: 


FWHM = K 


2cT 


( 8 ) 


where K is a constant, depending on the design. The constant, K, is always 
greater than one, and for most designs will be less than three. 

The finesse, ^ of a filter is the ratio of its free spectral range to 
its bandwidth. Since the FSR is determined by the greatest common divisor 
of its crystal lengths and the FWHM is determined by the total crystal length, 
the finesse is a function of their ratio. Specifically, 




2T 

KAt * 


C 9 ) 


For the Lyot and Sole filters. 


^Lyot ‘ 1 ' 13 - 2n 
^Solc = 1 ‘ 26 ‘ n 


CIO) 


where n is the number of elements in each, filter. 

To be comparable to a Fabry -Perot filter whose finesse is usually between 
20 and 50, a Lyot filter would have 4 to 6 elements and a Sole filter would 
have 20 to 50 plates. Other, more highly apodized designs, would have corres- 
pondingly more plates for a given finesse. Partial polarizing filters, for 
example, have about 25% more elements than a corresponding Lyot filter. 

Each of the birefringent filter types has the freedom to control trans- 
mission profile parameters such as the bandwidth and sidelobe level. Any 
change in the filter design which introduces a taper in the impulse response 
will have a corresponding effect on the transmission profile shape. There 
are, therefore, two major questions. The first is, what are the relation- 
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ships between desirable transmission profiles and impulse response 
tapers. The second is, what impulse response tapers are available with 
birefringent filters and how are they achieved. 

Although the first issue, the relationship of desirable spectral 
properties to impulse response, is an extremely sophisticated problem, 
it has received considerable attention recently. The impulse response , 
characteristics of birefringent filters are identical to those of finite 
impulse response digital filters and discrete array antenna systems. Conse- 
quently the results from these disciplines are directly applicable to the 
design of birefringent filters. An extremely complete catalog of known 
filter responses is given by Harris (1978) , though only selected designs 
and design criteria will be discussed here. 

The fundamental design tradeoff is between bandwidth and sidelobe 
(secondary maxima) level. In general filter designs which result in 
impulse responses which peak in the center and taper downward towards the 
edges result in lower sidelobes but wider bandwidths. This tends to 
concentrate most of the filter transmission in the main passband, but at 
the expense of some resolution. Two similar, but slightly different design 
objectives are, first, to minimize the peak sidelobe level for a given band- 
width, and, second, to minimize the integrated out of band transmission for 
a given bandwidth. The first criterion is natural when the filter is to be 
used for examining the line structure of emission spectra in the presence of 
very low background radiation. For examining the structure of absorption 
spectra in the presence of a large background continuum, the second criterion 
is more applicable. It has been demonstrated that the first objective is 
achieved by Dolph-Chebyshev tapers (Harris, 1978) while the second is 
achieved by discrete prolate spheroidal sequences (Slepian, 1978) . In 
addition, many useful suboptimal designs .are known which have other desirable 
properties . 

The question of which impulse response tapers are desirable may be 
considered to be solved for these purposes. This leayes the remaining 
question of which tapers are achievable By birefringent filters. In prin- 
ciple the answer to this question is that all tapers are achievable by bire- 
fringent filters. As discussed previously, Harris (1964), has presented a 
synthesis algorithm for achieying any imp.ulse response with a Sole filter. 
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In particular, Schiffman and Young (1968) have demonstrated that this 
technique can, in fact, be used to build Dolph— Chebyshev birefringent 
filters . Figure 1 shows the response of these designs for 15-plate 

Inn t 

filters with sidelobe levels of 10 ,10 , 10 and 10 . An interesting 

characteristic of these profiles is that, while the peak sidelobe level 
is minimized, all the sidelobes are equal. 

A corresponding synthesis for filters with perfect internal polarizers 
was presented by Ammann and Chang (1965) . It demonstrates that any speci- 
fied impulse response may be achieved by such a filter. Ihe disadvantage 
of this procedure is that the internal polarizers decrease the transmission 
substantially from that achievable with a lossless Sole filter. As an 
advantage, however, is the possibility of using fewer crystals, of differing 
lengths, such as in the Lyot design. The number of polarizers is not always 
prohibitive, therefore. 

Practical considerations impose some restrictions on the list of achiev- 
able responses. For Sole type, lossless filter^ the large number of elements 
required is the major problem. With larger numbers of plates, the angle 
tolerances for setting the plates become more severe. Present technology 
for producing precisely matched elements, coupled with the computational 
capability for initial alignment, makes these filters more practical now 
than in the past. Designs with more than 50 plates, howeyer, are probably 
only speculative at this time. 

Perfect internal polarizing filters are limited by the number of 
required polarizers, and by the difficulties of manufacturing elements of 
different but precisely matched lengths. As will be discussed subsequently, 
tunability allows for the compensation of length mismatches, so that problem 
may be avoided. The transmission loss of polarizers, however, presents a 
limitation in the use of these filters for arbitrary spectral transmission 
profiles. The solution to the practical limitations has been to use 
feasible design techniques which result in subop timal, but satisfactory, 
profiles. A very common design modification for Lyot filters is the intro- 
duction of a contrast element, an element whose length equals the second 
longest element in the original Lyot filter. This introduces an amplitude 
taper which reduces the sidelobe level significantly. In Figure 2 a four- 
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element Lyot response is compared with, a five-element filter which, uses 
a contrast element. Also shown in Figure 2 is the response of an 
untapered filter whose total length is equal to the contrast element 
filter. 

Partial polarizing filters, having more degrees of freedom, certainly 
have all the capabilities of both perfect polarizer and no polarizer filters. 
Their utility is that because of the extra degrees of freedom, more practical 
realizations for filters may be found. In particular, designs which use Lyot 
element lengths or nearly Lyot lengths but which haye better sidelobe control 
have been the major interest. The Lyot element lengths are the most efficient 
since they represent the division of a total crystal length into the fewest 
pieces while achieving a given minimum length, cAt. Any division using fewer 
plates will have zeros in the interior of its impulse response. 

The simplest partial polarizing design is one in which the first internal 
polarizer of a Lyot filter is replaced by a partial polarizer with: 

P = p y /p x = 1/3 • cm 

The response of such a design is very similar to that of a contrast element 
Lyot filter . A four— element Lyot filter with one partial polarizer is shown 
in Figure 3. The response of that filter should be compared with the contrast 
element Lyot shown in Figure 2. There is a 26% increase in bandwidth but this 
is due to the 27% more crystal in the contrast element filter. The sidelobe 
levels are comparable, however. 

If all the polarizers are replaced by partial polarizers in a Lyot 
filter the profile does not improve. There are, however, designs in which, 
more than one polarizer can be made partial and improvement is expected. 

These are discussed by Title C1976) as the alternate partial polarizing (APP) , 
double partial polarizing (DPP) , and redundant alternate partial polarizing 
(RAPP) , filters . 

The APP design is one in which, alternate internal polarizers are partial, 
p = .333, and the others are perfect, p = 0. Across a partial polarizer the 
crystal lengths are in the Lyot ratio, 2:1, but across a perfect polarizer 
the ratio 3:2 is needed. A four-plate alternate partial polarizer filter is 
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shown in Figure 4. As is typical of partial polarizer designs , the plates 
following the partial polarizers are orthogonal to the plates preceding the 
partial polarizers. Notice the extreme sidelobe reduction, while the band- 
width expansion is due to the reduced total crystal length. 

A double partial polarizing module may be designed with three Lyot 
length elements arranged in the order 02:4:1) with partial polarizers between 
elements . For this arrangement the partial polarizers should have orthogonal 
transmissions, p, equal to .333 and .111 respectively. Figure 5 shows the 
responses of a single DPP -module and, for comparison, the equivalent perfect 
polarizer Lyot filter. Notice that the sidelobe level has been reduced but 
at the expense of an increase in bandwidth. A filter can be constructed 
of double partial polarizing modules alternating with perfect polarizers 
just as an APP filter is constructed from single partial polarizing modules. 

A partial polarizing design which uses the idea of duplicating element 
lengths, like the contrast element design, is the redundant alternate partial 
polarizing configuration. The elements have Lyot lengths, but all but the 
longest and shortest elements are duplicated. The polarizers alternate, 
partial and perfect. Figure 6 shows a RAPP configuration with six elements, 
and the same configuration with perfect polarizers. Notice the excellent 
sidelobe reduction at only a modest bandwidth increase. In addition to 
their spectral profile advantages, the partial polarizing designs have been 
found to be much less sensitive to construction errors than the traditional 
Lyot, perfect polarizer, design. In particular, crystal length, errors produce 
much weaker secondary maxima in the partial polarizing designs than in the 
Lyot filter. Also, surprisingly, some partial polarizing designs have been 
found to be particularly tolerant to the polarization, p. This is particularly 
important for short wavelengths where sheet polarizing materials haye poor 
polarizing properties. Specifically the polarization of sheet polarizers 
deteriorates rapidly at wavelengths Below 4500 X so a design which is robust 
with respect to p is very desirable. Figure 7 shows the behavior of a partial 
polarizer design whose elements have lengths (4:8: 6:4: 2:1) . Notice that as 
p increases the finesse decreases, and the sidelobe structure changes, but 
the filter does not deteriorate catastrophically. A current topic is the 
search for partial polarizing designs which have their best performances when 
P is about .25. That should insure a greater uniformity over a broad range of 
p . Designs of this type should make it possible to design filters which are 
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acceptable down to 3500 2 using the same materials and construction techniques 
as used for longer wavelengths . 

The major transmission losses in birefringent filters are caused by- 
absorption by polarizers and imperfect filter construction. A minor source 
of loss is interface Fresnel reflection at the plates. Howeyer, for filters 
that are in an index matching fluid this problem is minimal. A ten percent 
transmission loss occurs for 100 calcite or 3400 quartz surfaces. 

All birefringent filters require high quality entrance and exit polar- 
izers . These can be ordinary sheet polarizer, but it is also practical to 
use thin film or crystal polarizers which can have virtually 100% efficiency. 
Unfortunately both thin film and crystal polarizers are nearly cubical and 
a serious length penalty is incurred if such, polarizers are used internally 
in the filter. The additional length introduces optical complications, so 
that for most designs sheet polarizing materials must be considered the only 
option for the majority of the internal polarizers. 

A perfect polarizer filter requires (n-1). internal polarizers for 
elements, so that relative to a lossless filter the transmission is 


n-1 

T = T 

Lyot o 


aii 


where is the transmission of the internal polarizers in the pass direction. 
For a partial polarizing filter 


n 

T Partial ~_J^ T i 

where the are the transmissions of the partial polarizers. Figure 8 
is a plot of total transmission versus number of internal polarizers for 
various r . This should be compared with Figure 9 which shows the finesse 
of a Lyot filter as a function of the number of internal polarizers . 

Figure 10 has plots of transmission, t , versus wavelength for high quality 
specially treated polarizer materials mounted in AR coated glass. Visible 
wavelength polarizers superior to those in "Figure 10 can be produced, but 

at present are not stable fop moye than a £ew months. Beyond 7500 2 there 
are at least two methods to produce good polarizers besides polarizing sheet. 
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WAVELENGTH (A) 


Figure 10 . Transmission versus wavelength of specially treated polarizing 
materials . 
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These are sheet crystal polarizers and wire grids. At present it would 
appear possible, using modern holographic techniques, to make wire grid 
polarizers useful for wavelengths beyond one micron. Below 4000 2, however, 
there is a serious problem with polarizers. UV polarizing sheet does exist, 
but its maximum transmission is in the region of 60%. It might be possible 
to create partial polarizer sheet material below 4000 2 with reasonable 
transmission, which would make the partial polarizer designs very attractive. 
Another possibility is the use of electron beam and X-ray lithographic tech- 
niques for production of wire grid polarizers. At the present state of the 
art it might be possible to create polarizers for wavelengths down to 1500 2. 

For the visible Figures 8 through 10 provide a useful guide to what can 
be done using Lyot type designs. If a total transmission criterion is at 
least ten percent transmission, then operation around 4000 £ is limited to 
a finesse of 64. At present the best blocking filters have FWHM of 6 2 in 
this spectral region, so that a minimum FWHM of the Sirefringent filter is 
0.1 2.. On the other hand it is possible to construct a filter with a finesse 
of 2048 which' operates from 4600 £ to 7500 & with a minimum transmission of 
35% and a maximum transmission of 64%. If the wavelength range is restricted 
to above 5000 2 ., a transmission of 48 to 60 percent can be achieved. 

At present it is reasonable to construct 0.1 X filters for the yisible 
and difficult to produce filters narrower than .05 2 . A 0.1 2 filter could 
be blocked for the entire range 4800 - 7800 2 . and haye a transmission in the 
range of 25% to 50% over the entire spectral band. For a more moderate resolu- 
tion of say 2 2 , the required finesse implies a transmission range of 35% to 
65%. And for a 25 2 filter, the transmission range would be 53% to 75%. 

Polarizer limitations make it necessary to consider other designs which 
try to combine the transmission efficiency of Sole filters with the crystal 
utilization of Lyot filters. As has been demonstrated in the previous section 
it is often possible to find Sole type filters with some multiple length plates. 
A particularly dramatic example of this is shown in figure 11. The nine-plate 
Sole filter with a triangular taper impulse response is shown, along with a 
five— plate compression. Notice that the first and last three plates have 
been replaced by single, length three, plates at their average orientations. 

The performances of these two filters are equivalent, but the compressed 
filter has only half as many plates. In Figure 11(c) this five-element 
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Figure 11. Transmission versus frequency for n 
approximation (±>) , two five-element 
and a six-element Lyot (d) . 
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ups with element ratios of 6:1 (c). 
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filter has been used as a module and duplicated using plates six times as 
long. The two sections are separated by a perfect polarizer which combines 
the large free spectral range of the short section with the narrow bandwidth 
of the longer section. The resulting filter compares favorably with the 
six-element Lyot filter which has the same total crystal length. It has 
ten elements instead of six, a finesse of 49. instead of 72., but it only 
has one internal polarizer instead of five, and sidelobes which are a 
factor of 50 below the Lyot. This or similar filter designs would be very 
attractive for the ultraviolet or infrared wavelengths where polarizers 
are critical. 
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(4) Tuning Birefringent Filters 

There are two methods for tuning birefringent filters, a direct method 
which is conceptually simple but difficult in practice and an indirect 
method which is practical but difficult to understand. The direct method 
is to actually change the optical length of each crystal plate in order 
to make each plate have a retardation which is a multiple of the wave- 
length of interest, while maintaining the proper length ratios. This 
may be accomplished by using pairs of wedge shaped plates as adjusters 
for each element. By moving the wedges in and out, the thickness in the 
optical path may be changed and the elements -may be tuned. The optical 
path length may also be changed by controlled temperature changes . A 
change in temperature will change both the physical length and the bire- 
fringence. Although direct adjustment is the method for tuning air spaced 
Fabry— Perot and Michelson interferometers it has not been the most practi- 
cal method for birefringent filters. 

In practice most birefringent filters are tuned By rotation of thin 
waveplates or polarizers. The key to understanding this method of tuning 
is to consider the effect of a birefringent plate on the polarization 
state of monochromatic light. Thus spectral or frequency domain analysis 
is more suited to this problem than the pulse response or time domain 
techniques of previous chapters. 

Consider linearly polarized light incident on a Birefringent element 
as in a Lyot filter. The plane of polarization is assumed to be at 45° 
to the fast axis of the plate and the incident light is monochromatic 
with wavelength, X. In the coordinate system of the plate the incident 
light may be expressed as : 


A, 


1 c 

— cos C2ir — t) 

/2 X 


,.i 1 c . 

A = — cos C2ir — t) . 

s /2 X 


a> 


The birefringent plate introduces a pure time delay. At, in the slow axis, 
which is equivalent to a wavelength dependent phase shift 
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A(j> = 2trc . 

A 


The light exiting the plate is, therefore, represented as 


A® - -z cos (2tt ~ t) 

/2 A 


A® = cos (2ir 7 t - A<j>) . 
s /2 A 


( 2 ) 


(3) 


When the phase shift, A<j>, is a multiple of 2ir the two components of 
(3) are in phase and the output is linearly polarized parallel to the 
incident plane of polarization. If A<j> is an odd multiple of it, the output 
is again linearly polarized, but in a plane perpendicular to the input 
plane . For other values of A<f> , the output light is , generally , ellipti— 
cally polarized. The axes of the ellipse are parallel and perpendicular 
to the input polarization plane. Rotating the output, coordinate system by 
45° to express the output in terms of components parallel and perpendicular 
to the entrance polarizer: 


A = A , + A 
|j f s 

A 1 =_A f + A s • 

After some algebra (3) and (.4) imply: 

A|j = cos ^ cos ( 2 n f- t - fi) 

= sin sin (2* ~ t - y^.) . 

Equations (5) may be combined, eliminating the time dependence, 
the ellipticity 


(4) 


(5) 


to exhibit 




( 6 ) 
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Equation (6) is, as asserted, the equation of an ellipse whose axes are 
parallel and perpendicular, respectively, to the entrance polarizer, but 
whose ellipticity is a function of the wavelength. 

In order to tune this element, therefore, the plate must be followed 
by an optical element which is sensitive to the ellipticity of the output. 
A polarizer is a detector of linearly polarized light, and is therefore 
sensitive to the limiting cases : 


- 2niT 


C7) 


A<J>£ = C2n-l).ir . 


Using a polarizer and switching from the parallel to the perpendicular 
orientation switches the wavelength from 


to 



. _ 2cAt 

A 1 2n-l 


( 8 ) 


A more general ellipticity sensitive device is cons true ted “by insert- 
ing an achromatic quarter waveplate between the element and the polarizer 
as shown in Figure 1. An achromatic quarter waveplate is a birefringent 
plate which behaves as implied by equations (3) but for which the phase 
shift, A<j> = ir/2, independent of wavelength. The theory of these waveplates 
is described in Title (1975) and McIntyre and Harris (1968). which, are 
included in Volume II, and in Pancharatnam (1955a and 1955b) and will 
not be described in this section. It is, however, the ease of making 
achromatic waveplates inexpensively which makes tunable birefringent 
filters practical. 

While the theory underlying the construction of achromatic quarter 
waveplates is difficult, their application as tuning elements is not too 
difficult. Consider linearly polarized light incident on a quarter waveplate 
at some arbitrary orientation, 0, with respect to the fast axis of the plate. 
In the coordinate system of the plate the incident light is represented as: 
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Figure 1. Diagram of simple tuning element. 
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= cos 0 cos C2ir t) 

( 9 ) 

A"*" = sin 0 cos (2ir ~ t) . 

s A 

The light exiting the quarter waveplate is, 

A® = cos 0 cos C2ir ^ t) 

A e = sin 0 cos C2 tt 7 t - {) (10) 

s A l 

6 C 

or A = sin 0 sin (2tt — t) . 

S A 

Equations ( 10 ) represent elliptically polarized light whose axes are the 
fast and slow axes of the waveplate. The ellipse so defined is represented: 



A quarter waveplate, therefore, converts linearly polarized light into 
elliptically polarized light. By reciprocity, however, it must also conyert 
elliptically polarized light into linearly polarized light. The crucial- 
point is that the plane of polarization of the linearly polarized light 
will depend only on the ellipticity as determined by (11) . 

The mechanism is, therefore, that each wavelength (.modulo a free 
spectral range) will emerge from the birefringent element as elliptically 
polarized light with axes parallel and perpendicular to the input plane of 
polarization. An achromatic quarter waveplate with its fast axis parallel 
to the entrance polarizer, at 45° to the fast axis of the element, will 
convert the elliptically polarized light to linearly polarized light at an 
orientation which depends on wavelength. Specifically, from (. 2 ) , (fj) , and 
(11) , the dependence is : 


cos0 = + cos 


2 


( 12 ) 
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or 



ircAt 

X 


ror . 


03) 


At 0 = 0 the natural wavelength, X q , of the element emerges, so (13) may 
be written 


9 - TCit (x - r ) «*'> 

o 

The wavelength, X, is then selected by rotating an exit polarizer to the 
appropriate angle, 0. 

The tuning scheme described above is the most commonly used system 
for Lyot filters. It has, however, one serious liability. Since the exit 
polarizer for the first element is the entrance polarizer for the second 
element, the second and subsequent birefringent plates must rotate as the 
filter is tuned. While rotating the crystals is theoretically satisfactory, 
as a practical matter it is very bad. No crystal elements are perfect, but 
the overall imperfections of a fixed set of crystals can be calibrated or 
polished out. If the crystals rotate, differentially, the errors become 
wavelength dependent and are yery much more difficult to compensate. 
Additionally, if the rotation axes differ from the optic axis, displacements 
occur which are also wavelength dependent. 

The problem of rotation of successive elements can be solyed by the 
addition of an achromatic half waveplate in the tuning element. An 
achromatic half waveplate is a birefringent plate which behaves as implied 
by equations (3). but for which the phase shift, A<i> = it, independent of 
wavelength. Conceptually they are just two quarter wave plates in series. 
When linearly polarized light is incident on a half waveplate at an 
orientation, 0, with respect to its fast axis it exits linearly polarized 
but at an orientation, 20. The half waveplate may, therefore, be used to 
rotate the plane of polarization in between the exit polarizer and the 
succeeding crystal. This configuration is shown in Figure 2. 

The only problem with the position of the half waveplate in Figure 2 
is that as the polarizer rotates to tune the element the half waveplate 
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EXIT POLARIZATION 
DIRECTION (0°) 


Figure 2. Diagram of a tuning element with a half wave plate behind it to 

rotate the exit polarization direction parallel to input direction. 
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must rotate half the angle In order to correct the orientation* If, 
however, the half waveplate is placed before the exit polarizer, as 
in Figure 3, only the half waveplate need rotate* The exit polarizer 
and the succeeding crystal will stay fixed. Each of these placements 
of the half waveplate have been used in the Lockheed tunable filters* 

The tuning technique for partial polarizing and Sole filters is 
more complex than that for Lyot filters. In Lyot filters the tuning 
only has to be concerned with one plane of linear polarization since the 
entrance and exit polarizers reject the other plane. In the context of 
the configuration shown in Figure 3, the element is tuned to a wavelength, 

X , if vertically polarized light at the entrance polarizer emerges verti- 
cally polarized at the exit polarizer. This is the behavior which xras 
demonstrated above. The Sole and partial polarizing filters, however, 
place more demands on the tuning section. They require linearly 
polarized light of wavelength, X, at any orientation, to emerge from 
the tuning section at the same orientation as the input. The polarizers 
are, of course, absent. 

In order to determine the effect on arbitrary orientations it is 
only necessary to examine the effect on horizontal and vertical polariza- 
tions . All other polarizations are simply linear combinations of these 
two orthogonal polarizations. By extending the calculations of equations 
(5) to include the effects of the quarter wave plate and the rotating half 
waveplate of Figure 3 the responses to vertically and horizontally polarized 
light may be determined. For a vertically polarized input, the yertical 
and horizontal components of the output are, respectively: 


a / 3*11" c f A(|) \ 

Ay - cos fey- t + 2*4 

Ah-° . 


as) 


Similarly, for a horizontally polarized input, the output components are 


Ay * 0 

A e - cos (^f t + f - f*) . (16) 
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EXIT POLARIZER (0°) 



Figure 3. Diagram of a tuning element with a rotating half wave plate tuner. 
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Equations (15) and (16) assume that the system is tuned to the wave- 
length, X, by rotating the half waveplate by an angle 

(17 

An examination of (15) and (16) indicates that vertically or horizon- 
tally polarized input becomes vertically or horizontally polarized output 
respectively. This is what is required, except for the 3 tt/ 2 phase lag 
indicated by (17) . This phase lag implies that arbitrary linear input 
polarizations will become elliptically polarized outputs . The addition 
of an additional quarter waveplate, oriented as shown in Figure 4, will 
solve this problem by introducing an extra ir/2 phase lag in the horizontal 
component. The resulting configuration becomes a tunable element which 
may be used in Sole or partial polarizer filters. With the addition of 
entrance and exit polarizers it will still work in Lyot filters, although 
the extra quarter waveplate is then unnecessary. 

There are filter configurations for which, modifications must be made 
to the above tuning system because of the introduction of wayeplates in 
the filter for other purposes. These usually only inyolve a fixed rotation 
or the addition of achromatic wayeplates. The principles described aboye, 
however, remain the same. 
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QUARTER WAVE 
PLATE (0°) 



Figure 4. Diagram of wave plate arrangement for tuning a Sole or partial 
polarizer element. 
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(5) Field of View 

One of the most significant advantages of the birefringent filter 
is its large field of view compared to Michelson and Fabry-Perot inter- 
ferometric filters. This advantage, which is between 50 and 300 means 
that a birefringent filter can gather 50 to 300 times more light than 
a comparable diameter interference filter. In order to achieve this 
field of view advantage, the elements of the birefringent filter must 
be constructed as widefield elements. The widefield construction was 
discovered by Lyot and was subsequently described by Evans. Recently, 
it was discovered that special -materials and alignment perturbations 
can increase the field by another factor of 3 to 10. A complete analysis 
of these field of view effects appears in Title and Rosenberg (1979) 
which is included as an appendix to this report. In this section a 
brief analysis will be presented accompanied by a summary of some of 
the important results . 

A widefield birefringent element is made by dividing an ordinary 
element in half, rotating the second half around the optical path by 90° 
with respect to the first half, and then separating the halves by an 
achromatic half waveplate. The half waveplate is oriented parallel to 
the entrance and exit polarizers , at an angle of 45° with respect to the 
first crystal. As discussed in the preceding section, a half waveplate 
rotates a plane of polarization by twice the angle between the input 
plane and the axis of the wayeplate. An impulse exiting the first 
crystal along the fast axis, at 0°, will be rotated by the waveplate 
to an orientation of 90°. It will, therefore, also traverse the fast 
axis of the second crystal. Similarly, a pulse along the slow axis of 
the first crystal will also traverse the slow axis of the second crystal. 
For a normally incident light ray the widefield element will, therefore, 
have the same retardation as a simple, narrow field, element. 

The reason that simple elements aren't wide field is that, for 
non-normally incident rays, the retardation of a plate varies with the 
angle of incidence, i, as strongly as in a Fabry-Perot. Unlike a Fabry- 
Perot, however, it also varies with the azimuthal orientation, 0, of 
the ray . In particular for rays at 0° azimuth, with respect to the 
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crystal's optic axis direction, the retardation decreases with increasing 
incidence angle, i as: 

E(i,0°) (!) 

o 

At an azimuth of 90° the retardation increases with incidence angle as: 


R(i,90 ) = R 


1 + 


. 2.1 

sin l 

2n n 
o eJ 


( 2 ) 


where n Q and n^ are the ordinary and extraordinary indices of refraction, 
respectively, and R q is the nominal retardation of the plate. For arbi- 
trary azimuths the retardation is given, as a second order approximation, 
as : 


R 


= R 1 - 
o 


sin 2 i 

2n 2 


2 

(.cos 0 


sin 2 0). 


with p 



(3) 


The azimuthal behavior is the key to making a widefield element. In 
the quadrant containing the optic axis, 0=0°, the retardation variation 
is negative with incidence angle, as in (1), while in the next quadrant 
it is positive, as in (2). Thus, in the widefield element with the 
crystals orthogonal, these positive and negative quadrants oyer lap each 
other. An off-axis ray at 0° azimuth, as it goes through the first crystal 
will have less retardation than a normally incident ray. As it continues 
through the second crystal it will have a relative azimuth of 90°, since 
the crystal is rotated and therefore have more retardation than a normally 
incident ray. Overall, therefore, the retardation variation will be 
reduced as these positive and negative quadrants tend to cancel. To 
find out what the residual variation is equation (2) can be used twice, 
once for each crystal. The result is 
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W 1 ’®) 


R 

o 

2 


-b 



which simplifies to 


- * 0 



[* 


sin 1 , 2„ . 2_v 

— ~ — (cos 0 - p sm 0) 

2n 

o 


1 Ceos 2 (0+90) - p sin 2 (0+90)) 
2n 

o 


• 2 . 
sm x 

4n 2 


. 2 . 
sm x 


4n 



(4) 


(5) 


There is still some retardation variation with, i, but it is considerably 

reduced from (3) . The factor (1-p) is very small for all known birefringent 

materials, always less than 0.2. For example, for calcite it is .11 and 

for quartz it is .006. Thus, for calcite there is about a factor of 20 

2 

difference between the coefficients of sin i in (3) and (5) . For even 
smaller index differences, the angular retardation variation - decreases 
even further, but for very small differences the second order approxima- 
tions (3) and (5), are no longer adequate. 

The actual expression for the single crystal retardation is: 


R(i,e) = | 




(6) 


in which 


a = sin i sin 0 
b — sin i cos 0 
d = crystal length 
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and the nominal retardation is 


Though this is difficult to work with, increasingly higher order approxi- 
mations have been developed . These become important for small index 
differences or relatively large bandwidth elements. 

The relative variation in transmitted wavelength is the same as 
the relative variation in retardation, 


SX 

X 

o 



(7) 


When the filter is used in an optimal optical system the field of view 
is set such that the wavelength shift of the extreme rays is equal to 
the FWHM at normal incidence. Thus, the field of view is limited by 



max 


FWHM 

X 

o 


( 8 ) 


Equation (5) is applicable whenever 

FWHM ^ P e~ n o ^ 

X _ 2 

o 2n 

e 

For smaller index differences or longer bandwidths the higher order terms 
must be used. This is illustrated in Figure 1 in which the incident angle 
limits are shown for quartz and calcite, as a function of the bandwidth 
or relative wavelength tolerance. In Figure 1, the bandwidth shown is 
the equivalent bandwidth when the element is used in a telecentric optical 
system since this is the usual application. A telecentric system insures 
that the rays have the same angular distribution at all points within the 


(9) 
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Figure 1. Incident angle and sine incident angle versus bandpass and relative wavelength 

shift in a cone of that incident angle (top) for quartz and calcite. Also shown 
is the incident angle ratio and throughput ratio with respect to F-P and MI versus 
bandpass and relative wavelength shift (bottom)* On both halves of the figure the 
2nd and 4th order approximation are shown dashed * N 
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field of view, but at the expense of increasing the filter bandwidth 
by a factor of /l~ over collimated, normally incident rays. 

It is well known that for air spaced Michelson or Fabry-Perot 
filters the relative wavelength tolerance is 


61 _ _ sin^ i 

1 ~ 2 
o 


( 10 ) 


If (10) is compared with (5) it is seen that, for equal wavelength 
tolerances, there is a solid angle, or throughput, advantage for bire- 
f ringent filters of : 


. 2 . 

sxn i 


BR 


. 2 . 

sxn x. 


2n 
o 

1-p 


FP 


Cll) 


This throughput advantage, and the corresponding incident angle ratio, 
is also shown in Figure 1 as a function of the telecentric bandwidth. 

Besides the throughput per unit aperture advantage of birefringent 
filters, there is also a significant potential advantage in aperture 
size. Currently calcite is available in crystals of 40 mm to 50 mm 
diameter, but quartz can be obtained with diameters up to 300 mm. Even 
more exciting is the possibility of building filters using oriented plastic 
film as a birefringent material. These could have apertures as large as 
1 meter diameter. 

Table I summarizes some of the properties of birefringent elements. 

The 'index difference, throughput advantage, and FWHM-length product are 
shown for various materials. Two plastics are included, PVA, polyvinyl 
alcohol, and PET, poly (ethylene terephthalate) . PET is more commonly 
known by the duPont tradename. Mylar. As the table indicates, the bire- 
fringence of plastic film is not unique and depends on the degree to which 
the material has been oriented. 
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Table I 


Material 

(n -n ) 
o e 

Throughput 

Advantage* 

FWHM- 

length** 

Calcite 

.1720 

47.5 

.727 

KDP 

.0454 

151.6 

2.76 

MgF 

.0118 

450.5 

10.58 

Quartz 

.0091 

813. 

13.7 

PVA 

.003-. 013 

370.-1500. 

10.-42. 

PET 

.05-. 20 

25.-98. 

2.5-6. 


* 


Throughput advantage = 


2n 2 n 
o e 


xi -n 
e o 


and is relative to an air spaced 


Michelson or Fabry-Perot filter. 


** FV'THM— length is in units of 


cm 


-12 2 

(10 m ) , and is measured at 


X = 5000 
o 
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(6) Bandpass Control 

In a preceding chapter the question of profile management and 
bandwidth determination was discussed, but all the alternatives offered 
there were design decisions before the filter was constructed* None, 
except perhaps rotating all the crystal elements of a Sole type filter, 
had any capability to adjust the bandwidth after the filter was con- 
structed* There are, however, applications for which the ability to 
control the bandwidth of a filter is useful* As the bandwidth is 
increased the transmission increases, so if there isn f t a need for 
particularly fine spectral resolution a wider bandwidth could be 
desirable. A design for which the bandwidth may be adjusted in use 
would, therefore, be advantageous. 

For a Lyot filter, the bandwidth can be adjusted by a physical or 
optical removal of the longest element. That would double the band- 
width without changing the free spectral range. Physical removal, 
while conceptually simple, has many practical difficulties and is 
considered an unacceptable approach. Optical removal is possible and 
the rest of this section will discuss several, both simple and sophis- 
ticated, methods which are feasible. 

If the filter is a narrow field filter, then each element is made 
from a single crystal (with, perhaps, a tuning section on its end). 

Optical removal of such an element is enabled by using two crystals, 
each half as long as the original. Ordinarily the two halves would be 
aligned with their fast axes parallel. Their retardations will, of course, 
add and the element will be active in the filter. To remove the element 
one half is rotated by 90°. The retardations will subtract and the 
element will effectively be removed. This has, however, two disadvantages. 
The first is that total removal is wasteful. It would be more useful if 
only part of the retardation is removed, thus enabling the element to be 
used as a contrast element for profile control in the wider bandwidth 
mode. This problem could be remedied by splitting the element into non- 
equal pieces, which would make the element effectively short but non-zero 
in the rotated configuration. The second problem is that, as discussed 
previously, rotating a large retardation crystal element is undesirable. 
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Fortunately, as with a tuning section, achromatic quarterwave and half 

waveplates may be used to accomplish the rotation. If a half waveplate, 
o 

oriented at 0 , is placed between two crystals which are oriented at 

+45 and -45 , respectively, then the retardations will add. Now, 

replace the half waveplate by two quarter waveplates, each oriented at 
o 

0 and nothing changes since they add to equal a half waveplate. If, 
now, one of the quarter waveplates is rotated to 90°, the quarter wave- 
plates cancel, and effectively disappear from the optical path. The two 
crystals, being orthogonal, also subtract retardation. By the rotation 
of a single quarter waveplate, therefore, the retardation of an element 
may be switched from long to short. Similarly, if the crystal is divided 
into several pieces, say of lengths 1, 1, 2, 4, 8, ... , separated by 
quarter waveplate pairs, then the element can be adjusted to synthesize any 
even length by appropriate waveplate rotations . 

In practice, however, the situation is almost as simple, but the 
theory must be extended to the wide field elements described in the 
previous section. Remember that a widefield element is constructed by 
splitting an ordinary element and placing a half waveplate between the two 
orthogonally oriented halves. The obvious realization of a variable band- 
width widefield element is, therefore, just two variable bandwidth ordinary 
elements, orthogonally oriented, separated by a half waveplate. As long 
as the two halves are each adjusted to the same retardation the combination 
will behave as a widefield element of twice the retardation. Thus, to 
make a widefield element which can be switched from a retardation of 16 to 
a retardation of 4 the following configuration would suffice: 

5(0°) :Q(45°) :Q(±45°) :3(0°) :H(45°) :3(90°) :Q(+45°) :Q(-45°) :5(90°) (1) 


where Q and H denote quarter and half waveplates, respectively, and the 
orientations of each plate are in parentheses. Such an element could be 
switched between a narrow bandwidth or a contrast element role in a Lyot 
filter 

The configuration shown above may be simplified, slightly by replacing 
each quarter waveplate pair with a single rotatable half waveplate. This 
doesn't work for a single narrow field element since the extra half wave 
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of retardation will shift the transmission, but in the split widefield 
element the two outer half waveplates will cancel. The resulting con- 
figuration will then be: 

5(0°) :H( 4 °o) :3(0°) :H(45°) :3(90°) :HC 4 °o) :5(90°) . (2) 

An important feature of this configuration is that its field of view 
properties correspond to its effective retardation, not its total crystal 
length. When it is adjusted to have retardation 16, it has the field of 
view of a widefield element of retardation 16 . After rotation of the two 
outer half waveplates, however, it has both the retardation and the wider 
field of view of a widefield element of retardation 4. Though this 
property is not obvious it will become more apparent in the subsequent 
discussion. 

Besides the variable bandwidth feature there is another reason to 
consider configurations such as the example above. In order to make the 
shortest elements of a Lyot filter, very thin crystal plates must be 
fabricated. Especially in widefield elements, where each plate is only 
half its normal thickness, this becomes a very difficult procedure. 
Additionally, the resulting plates are very fragile and impractical to 
use in any system which is subject to mechanical shock. For this reason 
it is desirable to use thicker crystals in the subtractive mode described 
above. For example, a widefield element of retardation 2 could be con- 
structed from longer crystals as: 

5(.0°) :4(.90°) :H,(45°) :4(.0°) :5(.9Q°) . (3) 

The remarkable property of configuration (3) is that it has the full field 
of view of a widefield element of retardation 2. It is, therefore, a true 
equivalent of such an element. 

There is a moderately simple extension of the mathematical analysis 
of widefield elements which allows for the analysis of more complicated 
crystal/waveplate configurations. In particular it specifies easily veri- 
fiable conditions for which a configuration is widefield. From the discus- 
sion in the Field of View section, the off-axis behavior of a single crystal 
is : 
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* * f /\ 

R(i,0) = R {1 - [1 - (1 + p)sin 0] } (4) 

° 2n 

o 

where R(i,0) is the effective retardation at the angle of incidence, i, 
and azimuth, 0, and p is the ratio of the ordinary to the extraordinary 
indices of refraction: 

P " V n e * C5) 

For the widefield element: 

R o C0°):H(45 O ):R o (90°) (6) 

the retardation of the two crystals add and the field behavior becomes: 


R^ = R(i,0) + R(i,0 + 90°) (7) 

2. 

= R {2 - [2 - (l+p)sin 2 0 - (l+p)cos 2 0]} 

2n 

o 

, 2 - 

“ 2R 0 " t 2R 0 - R o (l+p)sin 2 0 - R o (l+p)cos 2 0] . (8) 

o 

2 2 

Since the coefficients of the sin 0 and cos 0 terms are equal, the azimuthal 
dependence vanishes and (8) becomes: 




- 2R - R 


sln 2 i 

0 2n 2 
0 


[1-p] 


(9) 


The vanishing of the azimuthal dependence is equivalent to the widening of 

the field since p is nearly unity. In the following, the coefficients of 
2 2 

(l+p)sin 0 and (l*lp)cos 0 will be denoted by S and C respectively. When S 
equals C the element will be widefield. 

Suppose that the first half of the widefield element is split into 
two parts of lengths aR^ and (l-a)R o , and that the second part is oriented 


- 89 - 



LMSC/D674593 


at 90° ♦ The retardations of these two pieces subtract and the field 
behavior of the whole element becomes : 


R = aR(i,0) - (l-a)R(i,0 +90) + R(i,0 + 90) 


( 10 ) 


= a[R(i,0) + R(i,0 + 90)] 

. 2 . 

= 2aR o - aR Q — f - (1-p) . 
2n o 


(ID 


This is still a widefield element, but with retardation, 2aR . The field 

o 

behavior is determined by the coefficient of 


. 2 . 

sm i 


2n 


(1-P) 


and is easily seen to be that of an ordinary widefield element of retarda- 
tion 2aR . 

o 

It may be easily verified that for any division of the two crystals 
the field behavior is given by an expression similar to (11) . All that is 
required is that the sum of the crystal lengths on one side of the half 
waveplate be equal to the sum on the other side, assuming that all the 
crystals are oriented at ±45° with respect to the waveplate. If this is 
satisfied, then the element will be widefield. 

It is relatively easy to generalize this result to configurations with 
more than one half waveplate. Assuming that all the waveplates are at 45°, 
and that the crystals are at 0° or 90° the following rules characterize 
the behavior of a configuration: 

(1) Crystal at 0° 

If the crystal follows an even number of half waveplates it 
contributes positively to the retardation . 

If the crystal follows an odd number of waveplates it contributes 
negatively to the retardation . 
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The azimuthal contribution is to the coefficient, S, and is 
equal in magnitude and sign to its retardation contribution . 

(2) Crystal at 90° 

If the crystal follows an even number of waveplates its 
retardation contribution is negative . 

If the crystal follows an odd number of waveplates its 
retardation contribution is positive . 

The azimuthal contribution is to the coefficient, C, and is 
equal in magnitude and sign to its retardation contribution. 

If, after applying these rules to all the crystals in the element, the 
coefficients, S and C, are equal the element is widefield. Moreover, the 
field characteristics will be that of an equivalent widefield element of 
the same effective retardation. 

As an example of the application of these rules , consider the 
configuration: 

1^(0°) :E(45°) :L 2 (0°) :H(45°) :L 3 (90°) :H(45°) :L 4 C90°) . (12) 

The contributions to the retardation and the azimuthal coefficients are, 
respectively: 


R = L 1 " L 2 " L 3 + L 4 


S = 


L l- L 2 


C = 


~ L 3 + L 4 


The widefield condition is, therefore, 


(13) 


L 1 L 2 “ L 4 ' L 3 

and the retardation is, in the widefield mode. 


(14) 


2(L X " b 2 ). 


(15) 
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The two outermost waveplates in (12) may be optically removed by rotating 
them to 0°, as in (3), which results in: 

1^(0°) :L 2 (0°) :H(45°) :L 3 (90°) :L 4 (90°) . (16) 

The respective contributions are: 

R = L x + L 2 + L 3 + L 4 

S * L x + L 2 x (17) 

C = l 3 + l 4 . 

For this configuration, the widefield condition is: 

Lj_ + L 2 * L 3 + L 4 . (18) 

In order to use this element as an adjustable element, both (14) and (18) 
must be satisfied simultaneously. These imply 


L 1’ L 4 
l 2 - L 3 

with the retardation switchable between 


(19) 


R = 2(L X - L 2 ) 
and R = 2(L 1 + L £ ) . 


( 20 ) 


This is the same result as in (3) , but now it should be more apparent why 
the field behavior is the same as the equivalent widefield elements. 
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(7) Engineering Considerations 

The previous chapters have been largely theoretical and have examined 
the optical principles of birefringent filter systems. This theoretical 
development , however, has largely been guided by a desire to design opera- 
tional filter systems. It has been the goal of the entire filter study 
program to construct filters which are suitable for operational ground or 
space telescopes. Consequently, the concepts of the previous chapters 
are the theoretical foundations for systems which actually work. 

There are some engineering problems which are encountered in the 
reduction to practice of birefringent filters. Some of these are typical 
of optical systems, and others are peculiar to birefringent filters. 
Generally the engineering difficulties may be partitioned into three 
categories: mechanical, thermal, and electronic. While none of these 

present serious engineering obstacles to filter construction and operation, 
they represent legitimate considerations . 

The primary mechanical problem is mounting the crystals sufficiently 
firmly for the optical requirements, subject to the constraint that they 
are not fractured by shock or acoustic loading. A method of solving the 
mechanical mounting problem is described in detail in Appendix 2, Devel- 
opment of Filters for Spaceflight. The basic concept is the mounting of 
the crystals in a primary holder with a thin non-hardening adhesive such 
as silicon rubber. These primary holders are then mounted in a fixed 
(non-rotating) structure which serves to provide both mechanical alignment 
and good thermal contact to the crystals. As stated previously, both 
spectral and image motion problems occur when crystals are rotated, so 
for high quality imaging systems the crystals are fixed, and only thin, 
quarter and half wave plates are rotated. 

Before the development of inexpensive electronic computers, tuning 
narrow band filters presented a severe mechanical problem. All tuning 
elements of the filter were geared together and rotated by a single drive. 
Since these were Lyot filters with 2:1 crystal length ratios, each element 
had to rotate twice as far as its predecessor. Additionally, the optical 
lengths of the crystals had to be very precisely 2:1 or the system would 
quickly lose spectral alignment during tuning. These restrictions are 
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so severe that mechanically geared, narrow bandwidth filters never 
operated over more than ±20 £ without some differential corrections 
to the element cascade. 

With presently available minicomputers, however, these problems 
are easily solved. It is straightforward to tabulate, or have an 
accurate analytic expression for, all the transmission peaks in the 
channel spectrum of each filter element. Each element may then be 
tuned individually by rotation of its tuning waveplate , under computer 
control. To tune an element through its free spectral range requires 
a 180° rotation, implying a wavelength accuracy of FSR/180° per degree 
of rotational accuracy. Since the free spectral range of each element 
is just twice its full width at half maximum, a wavelength accuracy of 
FWHM/100 implies a rotational accuracy requirement of only 0.9° for 
even the thickest element. This is well within the capabilities of 
even modest shaft encoders . 

The rotational accuracy requirements for birefringent filters should 
be compared with the equivalent requirements for Fabry-Perot or Michelson 
interference filters. There, the path lengths must be adjusted to accura- 
cies of small fractions of wavelengths. The rotational requirements of 
tunable birefringent filters represent very simple mechanical problems by 
comparison. 

Thermal control of birefringent filters has been a serious problem 
in the past, but, like tuning, it is solvable with minicomputer control. 
The temperature sensitivity of a 0.1 X FWHM calcite element is about 
0.3 S/°C which means that holding the element to an accuracy of FWHM/10 
implies a temperature control to ,03°C. While this has been an extremely 
difficult requirement previously, it may be solved by using the filter's 
own tuning capability. Instead of trying to maintain a constant thermal 
environment to .03°C, it is easier to monitor the temperature and compen- 
sate for variations by- retuning each element. This requires precise 
knowledge of the temperature sensitivity of the filter materials , good 
temperature sensors, good thermal design of the filter housing and a 
small minicomputer to do the necessary calculations. Fortunately, none 
of these represent severe restrictions . The filter itself can be used 
to determine the temperature sensitivity of the materials by placing the 
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filter assembly in an oven and monitoring the spectral profile at various 
temperatures . The required sensors and minicomputer are now readily 
available., Two 0*1 X FWHM filters have been built using the temperature 
compensation by tuning concept and they have proven to be accurate to 
.01 £ over a temperature range of 40°C. 

The filter electronics package serves as an interface between the 
filter and the minicomputer. It must be able to operate the drive motors 
for the tuning sections and also monitor the temperature sensors. The 
electronics package is well within the state of the art and represents 
no serious problems. 

There are, of course, problems associated with the preparation and 
assembly of the filter components. These filters are sophisticated 
optical assemblies and must be given the same care as any precise optical 
device. Fortunately, as in other optical construction, there are highly 
sensitive interferometric testing techniques which can insure proper 
preparation and alignment. Moreover, each component or subassembly of 
a birefringent filter is, itself, an interferometer. Thus each component 
need only be viewed in monochromatic light to produce the diagnostic 
fringe patterns. 

The same analysis which was used to determine the field characteristics 
of a well— aligned wide— field element can also be used to determine the 
behavior of misaligned or mismatched elements. This has been done, and 
a catalog of fringe patterns for various states of adjustment has been 
prepared, and is included in Appendix 1. These representative patterns 
make it easy to identify and correct problems during the optical polishing 
and assembly stages. Each element must be prepared so that the retardation 
is uniform to within a fraction of a wavelength. This is not an especially 
difficult optical task since it is the retardation which must be uniform 
to X/10, not necessarily the crystal surface. Because the birefringence 
is much less than unity this makes the tolerance requirement weaker than 
the corresponding X/10 requirement for good optical .surf aces . In particular, 
it is weaker by a factor of about 10 for calcite or 200 for quartz. 

The distinction between optical surface quality and retardation 
accuracy is very important in birefringent filter construction. For a 
quartz element, an optical surface which is uniform to X/10 implies a 
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retardation which is uniform to X/2000. Ih Fabry-Perot filters, for 
comparison, there is no distinction. In order to achieve the required 
A/ 200 retardation accuracy it is necessary to actually polish the 
surfaces to that figure. Thus, a birefringent element can achieve more 
precise retardation uniformity with less stringent optical polishing 
requirements . 

Since it was recognized that there was no need to rotate the 
crystal elements themselves' it has been possible to completely separate 
the preparation and assembly difficulties from the operational considera- 
tions. It has been found that a properly aligned filter, once assembled, 
remains stable, to the extent of being space qualifiable. This is in 
contrast with many optical systems in which the critical operational 
parameters are the same as the critical alignment adjustments. For 
example, in a Fabry-Perot filter a critical factor is the plate spacing, 
but this is just the dimension which must be changed to tune the filter. 
In the birefringent filters the critical alignment parameters are the 
crystal element orientations, but. these are never changed during opera- 
tion. This insures the stability of these systems. 

Examples of narrow band tunable birefringent filters are shown in 
Figures 1 through 5. In Figure 1 all the optical components of the 
Lockheed Alternate Partial Polarizing Filter (LAPPU) are shown prior 
to assembly. The widefield elements, clearly showing the length differ- 
ences, are in the background. The rotatable half waye tuning plates, 
the perfect internal and external polarizers and the internal partial 

* Ji 

polarizers are in the foreground. An exploded view, indicating the 
assembly of these components is shown in Figure 2, and the final filter 
with tuning motors is pictured in Figure 3. The result of eyolution of 
the design of these filters is shown in Figure 4. This is a Spacelab 2 
flight prototype filter with more compact motors and a somewhat more 
sophisticated mechanical drive system. Fundamentally, however, it is 
very similar to the LAPPU filter. Each filter has a 0.1 1 bandwidth and 
an aperture of 35 mm. 

Figure 5 shows a view of a 'tunable 10 2 Lyot filter with rotating 
crystal elements. The aperture is 75 mm. Elements for this filter were 
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Figure 1. 


Optical components of the 0.1 X LAPPU showing wide field elements, 
half wave plates, polarizers and partial polarizers. 
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Figure 3. Complete LAPPU filter. 
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Figure 4. Flight prototype 0.1 X filter for Spacelab 2. 
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Figure 5. A 10 X tunable filter with 75 mm clear aperture. 
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manufactured using techniques normally used to produce silicon wafers 
for integrated circuits. Although the filter operates spectrally, 
there are some undesirable variations introduced by the rotating 
elements. The filter will be rebuilt to tune by rotation of wave- 
plates instead. It does demonstrate, however, that direct 2:1 gearing 
can be used for a filter of modest bandwidth. This filter is tuned by 
a single motor, although two motors can be used to increase the tuning 
speed between widely separated wavelengths. 
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(8) Conclusions , Projections and Recommendations 

Tit the start of this program it seemed clear that it would be 
possible to build filters, tunable through the visible spectrum, and 
with significant throughput gains over conventional interference filter 
techniques. That such filters could have bandwidths as small as .05 8 
or as large as desired also seemed evident. The 0.1 2 and 10. 8 tunable 
space flight filters which were built during this program have verified 
these beliefs. The theoretical and experimental results of the program, 
moreover, have indicated that even more may be expected. 

The detailed analysis of the field of view characteristics 
reported in Appendix 1 has shown that an additional factor of from 3 
to 10 in throughput is even possible. The field behavior predictions 
of this analysis have also been shown to be an effective diagnostic 
tool for the proper alignment of the widefield elements of the filters. 
Knowing the relationship between the field pattern and the alignment 
status makes it possible to determine the exact state of an element. 

This makes the actual construction of filters much more practical. 

The bandpass control techniques which are discussed in Chapter 6 
indicate that a single filter can be an extremely versatile -instrument . 
Its operating characteristics may be selected, in use, to match the 
immediate needs of an observing program. Because an element can be 
configured so that it always has the field characteristics of its band- 
width equivalent, angular throughput may be traded off against spectral 
resolution while still exercising control over the out of band spectral 
profile. Thus, instead of a filter having a single operating character- 
istic which is matched to one scientific task it may have a selectable 
set of characteristics making it suitable for many tasks with no perform- 
ance degradations for any one. Such versatility can be especially 
valuable for space platform applications. 

Outside of the visible band, above 7500 & or below 4500 & the poten- 
tial of birefringent filters wasn't as clear at the beginning of this 
program. The most immediate problem was the availability of suitable 
polarizers. In the infrared, up to about 4 pm, sheet polarizer is 
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available. Though similar in operation to visible wavelength polarizers 
its efficiency is too low to make Lyot filter designs practical. Another 
choice, however, now seems promising. Modern holographic and shadow 
deposition techniques make wire grid polarizers a real possibility down 
to below 1 pm and perhaps as low as 7500 £. Several wire grid polarizers 
for the 1 pm region have been built and the process seems promising. An 
advantage to wire grid polarizers is that they have the possibility of 
100% polarization efficiency. 

For the shorter wavelengths, below 4500 £, the prospect for a 
birefringent filter was not promising. Only crystal polarizers appeared 
practical because the transmission of sheet polarizers drops to below 
50% very quickly as the wavelength decreases . The roughly cubical shape 
of the crystal polarizers makes them very unattractive for a Lyot cascade. 
The added length to the optical path presents serious engineering problems 
if many of the polarizers are used. For this reason partial polarizer and 
lossless. Sole, configurations were considered. 

Although a 2800 £ Sole filter had been built by Fredga and Hogbom and 
operated successfully, there seemed to be many objections to that design. 

For example. Sole filters seemed difficult to tune with crystal wedges; 
they were difficult to manufacture because of the large number of plates; 
and they were difficult to align because of the sensitive, interacting plate 
orientations . 

It was finally realized, however, that the conventional waveplate' 
tuning techniques, discussed in Chapter 4, would apply to Sole and partial 
polarizing designs. This makes the Sole design seem very attractive. 
Similarly the silicon wafer manufacturing techniques were shown to apply 
to moderate aperture quartz elements by actually constructing a 10 S 
tunable filter from laminated quartz plates. Most recently, moreover, it 
has been realized that the synthesis algorithms developed by Harris ej: al. 
provide a solution to both the multitude of plates and the alignment 
problems . 

The surprising result of Harris, namely that there are many plate 
orientation combinations which generate the same filter profile, offered 
the possibility of multiple thickness plates in some configurations . As 
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was shown in Chapter 6, this actually occurs and results in nearly a 
factor of two reduction in the number of plates. Also, as discussed 
in Chapter 2, the Harris synthesis procedure provides an excellent 
diagnostic tool for the alignment of a Sole filter. Now it appears that 
Sole filters will offer a practical solution for the shorter wavelength 
region, and for verification, an actual construction project is just 
beginning. 

Since Sole type filters have the capability of almost arbitrary 
profile synthesis it is also possible to design hybrid, Solc-Polarizer- 
Solc, filters to achieve moderate to large finesses. The design pre- 
sented in Chapter 6, while not optimal, is certainly very promising. 

As these designs use only one or two internal polarizers, rather than 
six or seven, they are still practical and make tunable ultraviolet 
filters a real possibility. 

The filter development program has made significant advances in 
both the theory and practice of tunable birefringent filter design: 

(1) For the visible, 4500 £ to 7500 X, region extremely narrow 
band filters with transmissions exceeding 50% in polarized 
light have been shown to be practical and versatile. Two 
such filters have been qualified for space flight. 

(2) In the near infrared, wire grid polarizers have been built 
and appear quite practical. As they offer the possibility 
of 100% polarization efficiency it should be possible to 
build practical filters in the IR. Work is progressing 
slowly in this area but should be increased. The Sole- type 
filters are also being considered for the near infrared. 

C3) For the .ultraviolet region between 1500 2 and 4500 X there 
are two possibilities. Partial polarizing filters should 
work well down to 3500 £ with virtually the same techniques 
as were used for the longer visible wavelengths . New designs 
are being investigated but the principles of construction 
have been established. Below 3500 £, however, the recently 
discovered Sole-hybrid configurations seem most promising. 
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Currently a crystal set is being fabricated so that experimental 
filter construction can begin this year. 

In all, the prospects for continued development of practical, high-perform- 
ance filters for space platform observations from 1500 $ to 10 pm seem 
very promising. It is even likely that the infrared developments may be 
extended to the millimeter region for which birefringent Sole filters 
have already been constructed. 

In order to enable the continued improvement of widefield tunable 
birefringent filters the following development programs are recommended: 

(1) Holographic and electron beam fabrication of polarizers for 
the infrared and ultraviolet. 

(2) Artificial growth of large diameter calcite crystals. 

(3) New birefringent material investigation for the various 
spectral regions . 

(4) Electro-optical waveplate construction. 

(5) Theoretical extensions to the general synthesis procedures 
to include hybrid and partial polarizing configurations . 

and finally, 

(6) The support and encouragement of non-solar applications of 
birefringent filters to utilize their full scientific capabilities. 
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ABSTRACT 


The field of view characteristics of wide field birefringent (WFB) elements 
are compared to those of the Fabry-Perot (FP) , Michelson (MI) , and wide field 
Michelson (WFM) interferometers. Throughput gains of 50 to 300 or more with 
respect to the FP or MI are demonstrated. Further, it is shown that by proper 
choice of material, WFB elements can have angular characteristics identical to 
WFM interferometers. The properties of misaligned and mismatched half length 
WFB elements are calculated. It is shown that properly "Misadjusted" WFB elements 
can exhibit throughput gains with. respect to properly adjusted systems. Finally, 
a catalog of fringe patterns for WFB elements of different materials and differing 
angular misadjustment is presented. 
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I. INTRODUCTION 

Wide field birefringent elements offer significant advantages in field of 
view compared to Fabry-Perot and conventional Michelson interferometers. This 
advantage has been realized by solar astronomers since the invention of the 
wide field element by Lyot^ in 1933. However, the relationship between wide field 
birefringent elements and Fabry-Perot and Michelson interferometers has not been 
generally appreciated. Further, although the basic approach of Lyot^ and Evans ^ 
in describing the field characteristics of the wide field element was correct, 
the approximate expressions that they reported are only valid for a limited range 
of cases. In particular, they are not valid for crystals significantly less 
birefringent than calcite. 

In the past, birefringent filters have seen use largely as solar instruments. 
These devices have typically operated in a single line such as hydrogen alpha and 
had bandpasses in the range of 1 to 0.125A* For such systems the descriptions of . 
Lyot and Evans are adequate; however, even in these cases the properties of 
elements with misaligned or imperfectly matched element half lengths cannot be 
properly handled with extensions of the second order expansions. 

At present it appears that the significant gains in throughput, (fifty to 
several thousand) of birefringent elements, the easily implemen table 
techniques for tuning over 1.25 decades in 'frequency and their superior out of 
band rejection make tunable birefringent elements very attractive systems in 
applications which require imaging of extended sources or light gathering from 
diffuse sources. In the more conventional spectrum scanning applications bire- 
fringent systems do not require the intensive data processing of Michelson systems 
nor do they require that the spectral region under investigation be sampled 
uniformly. In addition birefringent elements are not sensitive to their orientation 
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with respect to gravity nor are precision motion controls required for scanning. 
Temperature sensitivity, which used to be a problem with high resolution systems 
(R>10 4 ), can now easily be controlled by temperature sensing and thermal 
compensation via the tuning elements 4 . 

Because of the apparent advantages of birefringent systems it is worthwhile 
to establish their relation to the more conventional Fabry-Perot and Michelson 
interferometers. This is done in Section II. In Sections III, IV and V the 
formalism for the angular sensitivity of a simple birefringent plate, a well- 
aligned wide field element, and misaligned and mismatched wide field elements are 
developed. Using the formalism developed in the earlier sections the throughput 
advantages of the wide field birefringent element are calculated in Section VI. 

The purpose of Section VII is to explain the fringe patterns produced by 
wide field elements of differing material properties and differing alignment 
states. Experimentally, it has been found that the principle areas of transmission 
loss in birefringent elements are due to wide field misadjustment problems. 

Section VII demonstrates that virtually all sources of misadjustment can be identi- 
fied, understood, and corrected by reference to tne fringe patterns. 
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RELATIONSHIPS BETWEEN WIDE FIELD BIREFRINGENT ELEMENT , FABRY- 
PEROT, MICHELSON, AND WIDE FIELD MICHELSON INTERFEROMETERS 

It is important to establish the connection between the wide field hire- 

fringent (WFB) element and the more traditional Fabry-Perot (F-P) , Michelson (MI) , 

and wide field Michelson (WFM) interferometers. This can be efficiently done in 

terms of the time delays introduced by the various systems, Shoxm in Figure 1 is 

an optical schematic of a simple birefringent element and a standard Michelson 

interferometer. 

The simple birefringent element consists of an entrance polarizer, a bire- 
fringent plate with its optic axis in the plane of the crystal face and an exit 
polarizer. The purpose of the entrance polarizer is to prepare a single state of 
polarization, since light originating from different polarization states do not 
interfere. Waves propagating through the plate parallel to the optic axis 
propagate at velocity c/n^ and those perpendicular at velocity c/n^, where 
n^ and are the ordinary and extra-ordinary indices of refraction. Since the 
entrance polarizer is at 45° to the optic axis, equal amounts of light are in the 
two components. After propagating along a length d of birefringent plate, the two 
polarized beams are separated in time by 


A = — (n -n ) 

O n p r\ y > 


CD 


o c x e o" 

where A 0 is the time delay for waves at normal incidence to the crystal. The exit 
polarizer being parallel to the entrance polarizer passes equal components of the 
two beams. Because the two beams arose from the same state and after passing 
through the exit polarizer they are indistinguishable, they will interfere. The 
intensity relation for two beam interference with time delay A Q is just 


I (X) = cos 


2 ( ^ 0 ^ , 

A 
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Figure 1. Optical schematic of a simple birefringent element and a normal 
Michelson interferometer. 
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where A is Che wavelength and c is the speed of light. 

Now consider the Michelson; the incoming beam strikes a 50/50 beam splitter 
and half the beam is reflected upward while the other half is transmitted. The 
reflected beam strikes a mirror, M, and returns to the beam splitter. The trans- 
mitted beam strikes M' and is similarly returned to the beam splitter where half 
the light is reflected in the same direction as the transmitted beam from M. If 
M and M' are not equal distances from the beam splitter, there will be a time 
delay 

2 , 

A o ” ^ - d M’) (3) 

between the two beams. Since the beams arose from the same initial states and 
after merging they are indistinguishable, they too will interfere. Again, the 
intensity relation is just that given by relation (2). 

I (A) = cos 2 

A 

Note that the beam splitter functions twice in the Michelson, first to divide the 

incident beam into two beams, and second to recombine the two beams from the 

/ 

mirrors. That is, the beam splitter serves functions identical to both the 
entrance and exit polarizers of the birefringent element. The difference in spacing 
of the two Michelson mirrors has the same function as the difference in the optical 
path lengths along the two polarization directions in the birefringent plate. 

There are, however, some angle sensitivity differences between the bire- 
fringent element and the Michelson. These arise from the off axis wave propagation 
characteristics of birefringent materials. For a Michelson the angular dependence 
of the time delay is given by^ 
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^] ,/2 

L tl Z -* » 


( 4 ) 


where n is the index of refraction of the Michelson optical path, and 


a = sin i cos 0 
b = sin i sin 0 , 


(5) 


where i is the angle the incident ray makes with the normal to M' and 9 is the 
angle the projection of the ray makes with the plane of the interferometer. The 
theta dependence has been written implicitly to illustrate the formal similarity 
with the birefringent plate. From Equations (4) and (5), however, the theta 
dependence drops out for the MI. Most MI operate in air, but the index of the 
optical path has been included because the Michelson angular dependence is the 
same as that of a Fabry-Perot or interference filter which is often filled with 
a dielectric material. Equation (4) can be rewritten as 



The equation is usually seen expanded to second order in sine i, or 
Since 


U 

A„ 


SA 


( 8 ) 


Where <$A and <5 A are the variations with angle of the time delay and wavelength, 
respectively, the relation between relative wavelength shift and incident angle 



, . 2 . 

1 sm i 

2 2 
n 


( 9 ) 
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i* ~ i/1^2 

A 0 ~ 2 V 


Equation (9) is the relation normally used for evaluating F-P, MI, or interference 
filter systems. Neglect of the proper expression (Equation (6)) seldom causes 


problems . 


For a birefringent element the angle dependence of the time delay® is 




2 h 2 I1/2 
-2+ r-2) 

o n e 


a 2q_b2"|l/2 > 


where the theta axis coincides with the optic axis of the birefringent plate. 
Equation (10) represents the time difference between an ellipsoidal (extraordinary) 
and spherical (ordinary) wave. Expanding Equation Q.0 ) to second order, a 
relation similar to Equation (7) can be obtained 


A - cos 2 6 - (1^2.) sin 2 0 | 


where A 0 is given by Equation (1) . There are several important differences 
between relations (7) and (11). First, when 


cos 2 0 — (— — )sin 2 0 
n 0 


n 1/2 

tan0 = + (_J=.) 

“ «0 


(13) 

n o A o 

that is, the field of view is infinite at particular angles of theta. More 
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importantly, the sign of the angle sensitivity changes between the sectors defined 
by Equation (12). Along the optic axis and the orthogonal axis the angular 
sensitivities are, respectively: 


and 



sin 2 i 

2 n c/ 


(14) 


61 _ sin 2 i 

* 0 q = 9 q 0 2n o n e 


(15) 


Along the optic axis the angular sensitivity is the same as a MI as is to 
be expected. Orthogonal to the axis, however, the difference in curvature 
between the ordinary and extraordinary waves reverses the angular sensitivity. 

The reversal of the sign of the angular sensitivity in orthogonal directions 
is the physical phenomena that allows construction of a system of much lower angular 
sensitivity. The WFB element (see Figure 2) is constructed by splitting a normal 
element in half and rotating the second half so that the optic axes of the two 
halves are orthogonal. In addition, a half wave plate is inserted between the 
halves with its axis along the bisector of the angle between the optical axes of 
the crystals. Since the effect of a half wave plate is to rotate a plane of 
polarization by twice the angle between the polarization direction and the axis 
of the wave plate, the effect of the wave plate is to rotate, for states of 
polarization, the optic axes of the two halves back to parallel. Hence, the on 
axis time delay of the WFB and a simple element are the same, but for off axis rays 
the time delay is 


a - d 

A = -7T- 

2c 


{ n < 


i - c 


+ fc{- 


i - ( 


v 

i 2 


♦&> 

n e 


1/2 


n. 


n e 


- n 


-n r 


i -<£& 


n r 


1/2 


1 /- g-' 2 + b l2 v 

1 ( ) 


J 


1 / 




(16) 
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Figure 2* Optical schematic of a wide field birefringent element and a 
wide field Michelson interferometer. 
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where 

a' = sin i cos (0 + 90) (17) 

and b' = sin i sin (0 + 90), 

To the second order 

A " A o Cl - t— jr ( n e ~ n 0 ) sin 2 i) (18) 

HI1 o e 

or 


6A 


1 

4n 0 Z 


(^ 



sin 2 i . 


(19) 


Since the difference between the ordinary and extraordinary indices is less than 
one, the wide field construction has the desired result of decreasing the angular 
sensitivity. 

Note that Equation (19) implies that by decreasing the index difference, 
the angular sensitivity can be indefinitely decreased. This is not the case 
because of the higher order terms in the expansion. It is, however, worthwhile to 
consider the WFM, an analogous system whose field effects are determined by- fourth 
order terms. The wide field Michelson, shown in Figure 2, achieves its larger 
field of view by a block of material placed in front of M', of length, d^, , such 
that the image of M 1 , through the b lock, is optically coincident with M. The 
time delay 7 is 


A = 


. 2 . 1/2 . 2 1/2 

n a (“\ sxn i, , . /. sm 

" d H 01 ~ - n d M' G 7F" 5 

n n 


• (20) 
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To the fourth order 

* = • ' ( 21 ) 

Equation (21) was carried to fourth order because the second order term 
is zero under the wide field condition which occurs when 


d M* 



( 22 ) 


so that 


A = 


A o 


(1 + 


* 4 • 
sin i> . 

8n 2 n 7T; ’ 


(23) 


where 

2 

A o = ~( nd M - n ' d M'), (24) 

and thus 


6A _ sin 4 i f25l 

A 0 " 8^^ - 


If the lengths of the WFM are not exactly adjusted, so that 


d M' 



d M + 



then 


6A _ sin^j sin 4 i 

A 0 2n n' E 8n 2 n ' 2 


(26) 


(27) 


Clearly for the case of the WFM the field of view does not become infinite 
as e goes to zero, but is limited by the fourth order term. A similar effect 
should be expected for the WFB. It will be demonstrated below that for a well 

aligned WFB the angle sensitivity is virtually identical to an appropriately 
misadjusted WFM. 
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Iix. FIELD EFFECTS OF 'A SIMPLE BIREFRINGENT PLATE 

In order to analyze the field effects of not only perfect WFB elements, 
but also elements which are misaligned or mismatched in length, it is necessary 
to examine the higher order properties of a simple birefringent plate. The 
expression for the time delay, equation (10) ? will be rewritten as: 


4 - V{H 1_X2<1 ' »-P 2 )2>] 1/2 - £ [l - 

in which 




(a) 

p = 

n o 

rig ^ 

(b) 

q = 

i-p » 

(c) 


sin i 


x = 

y 

n 

(d) 


o 


z = 

sin 2 0 • 

(e) 



Symbolically this "is written • 


A = A 0 f (yr, p, q, z) (30) 

and has the power series expansion 

A = A 0 1 + c x x 2 0(1) + C 2 x 4 0(2)+...+C n x 2n 0(n)+. .. (31) 

The coefficients, C n , are the coefficients of S n in the expansion of 
(1-S) 1 / 2 . Specifically 


2 (_2n-2) ! 

^n n ! (n-1) ! 


( 32 ) 
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th 


The factors, 0(n), are n degree polynomials in z. From the expansion of the 

1/2 

form (1-S) it may be seen that 


0(n) = 


(l-(l-p 2 )z) n - p 


The first few terms of .expansions C32) and (33), are given by: 


(33) 


Cj = -1/2 


c 2 " - 1 / 8 


.2 Z 


0(1) = l-(l+p)z 

0(2) = l-2(l+p)z+q(l+p)V 

c 2 ” ”1/16 0(3) = l-3(l+p;z+3q(l+p) 2 2 2 q 2 (l+p) 3 z 3 


(a) 

(b) 

(c) 


(34) 


Before going further it is, worth noting one point, namely that 0(2) is 


nonzero when 0(1) is zero. Clearly 0(1) is zero when 


z = 1 


1+p 


or when 


At this point 


n_ !/ 2 
tan 9 = + (— 1 “) 

~ v n 0 


(35) 


(36) 


0(2) = -p . 

Therefore, to fourth order, for tan 2 0 = (^ £ -) , 

n o 


(37) 


6X _ sin 1 *! 

~ 8^7^ P • 


(38) 


The fourth, order term indicates a deviation from the straight line asymptotes at 
2 

tan 0 = + n e /n o> There are still arbitrarily large angles of incidence for which 
the time delay sensitivity is zero, but they occur at different thetas. 
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XV. FIELD EFFECTS OF A WELL ALIGNED WIDE FIELD ELEMENT 

The analysis of the single plate may now be used to calculate the 
properties of a perfect WFB element. Using the formalism of (30), the time delay 
of a WFB element is 


A = f°- | f(x, p, q } Z) + f (x, P, q, z') 
where, as before. 


(39) 


z = sin 2 0 


and 


z = sin 2 (0+90°) = cos 2 0 


(40) 


(41) 


As before, expression (39) may be expanded in a power series in the 
variable x ; 

- 2n 

A = A q l+C^-x 0 ^t f (1)+C 2 x 4 0 wf ( 2)+. . .+C n x °WF^ 2n ^+-*- 


(42) 


The coefficients, C n , are the same as (32), the coefficients of S n in the 
expansion of (1-S) 1 / 2 . The factors, however, are more involved for the 


wide field element. They are: 

. r ] 

°WF (n) = J^ |U-(l-p 2 )z) n + (l-(l-p 2 )z') n -2pj . 

In this expansion of A, the angle of incidence dependency is represented by the 


(43) 


power series in x. The theta dependency, the variation of field width as a 
^function of azimuthal position, is included in the polynomial factors, 0^ F (n). 

As will be shown, each of these factors has a dominant term which is independent 
of theta, and smaller terms which are functions of sin 20. Specifically this 
theta independent term will be the sensitivity along the 0=0° ray, and is 
particularly easy to evaluate. 
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For this special case we have 


z = sin 


z 1 = cos^ 0=1, 


and (43) reduces to 


°WF (n) I G = 0 


- 2p + 1 . 

2q 


To sixth order in the sine of the angle of incidence, the sensitivity is given as: 

£■- T ^ - il - p 2 aw fP - i . (46) 

o n o n o 

Since q is small and p is a number near unity for most crystals, the factors, 

<V»> i - 

© = o l 1 - n. n > 1 


Along the ©= 0 ray, therefore, the angle sensitivity is approximately given by 

2 * / s 

§2l s - i s ^ n i + 1 sin i , 1 sin i (&&} 

A 4 2 8 4 + 8 6 * 

° “n ^ n 

o o o 


Equation (48) resembles, at least to fourth order, a slightly misadjusted WFM 

(equation 27) . An interesting consequence of equation (47) is that the coefficients 
n 

of P 1 decay slowly. In particular, 

C n C W ( ^ I ” (n-1) n>l (49) 

0=0 4 n n! (n-1) ! 

= 2 f 2n-2 | 

,n n J • 


The first few such factors are 




e= o 


= .125 

n = 2 

= .125 

n = 3 

= .117 

n = 4 
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The consequence of this is that for moderate field widths low order expansions 
will not he accurate* 

Xn order to consider the general form, (43 ) y including the theta dependency 
it is convenient to expand 0^(n) using the binomial theorem. This results in 

n 

°WF (n) = 1 + Iq 23 ( “ 1)1 (1-p2)1( i) (s^ 21 0 + cos 2l 0) n>l. (51) 

i=l 


The factors 


- . 2i_ . 2i 

= sm G + cos 6 


may be expressed as 

i/ 2 

S ± = 1 - / , A ± Ck) sin 2H 2G 


C52) 


(53) 


in which the coefficients, (k.) , may be determined from the relations: 

i>2 

A i+l (k) “ * k.>l, i>2 


\ (1) = i/4 


(54) 


The first few expressions for the sums, S , are 

i 

= sin^G + cos 2 ©' = 1 


S 2 = sin^G + cos^G - 1 - ~ sin 2 2G 
S, = sin^G + cos^G = 1 - y- sin 2 2G 


(55) 


S S o 1 A 

S 4 = sin G + cos 0=1- sin 2G -f* — sin 20 
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Substituting expressions such, as (55) into (51) enables the evaluation of the 

factors, Oyj,(n) , including the theta dependency. Again the first few such factors 
are: 


°WF^ 1 ^ = q/2 


o t ™( 2 ) = £--- - 2p+1 

WF V ' 2q 


( 1+p) 2 q sin 2 20 
4 


0 ttt ,( 3) p 6 - 2p+l 3 . ,2 2. 2 

W 2q " 8 a+p) q a+p > sin 20 ( 56) 

0 HF (4) . P 8 - 2p-KL - | (l+p) 2 q (l+p 2 +p 4 ) sin 2 2 a+ 1 U+p)V sin 4 20. 

2q -Lb 


For p very near unity, equations (56) become, respectively: 


° w <1) = q/2 


0 WF (2) - -1 -q Sin 2 2© 


°WF (3) » -2 -3q sin 2 2e 

0 W p(4) 55 -3 -6q sin 2 2G+ q 3 sin^G . 


(57) 


Each of the factors, O^Cn), may be seen to consist of a constant, theta independent, 
term and theta dependent terms which are even powers of sin2Q. This is the decom- 
position into the independent and dependent terms , for which the independent terms 
represent the behavior at 0=0 . Additionally, they represent the behavior at all 
multiples of 90° around the element. 

Although these theta independent components decrease slowly, as indicated by 

equations (50), the theta dependent terms decrease rapidly by comparison. This is 

because the first appearance of a sin 2n 20 term is in the 0yy(2n) factor. It is 
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therefore already multiplied hy x n , or sin* n i/n * n . In addition, the coefficient 
2 n 

of sin 20 in the expression for 0 TTO C2n) is: 

Wr 


C-l) n (l+p) 2n a ^ 
4 n 


158) 


For p close to unity this is approximately 

0-1) n< i 2n_1 

which decreases very quickly as n increases ♦ 


(59) 


The time delay expression, equation (42) , could then be written as 


A = i 

1 + A (i) 1 

+ A. (i, 0) \ 




0 1 

0=0 

/ ' 


(60) 


This implies that 

the sensitivity 

is decomposed as 

<$x sx 

X X 

I . AX 

1 0 = 0 + • 

(61) 


o o o . 


The first term on the right hand side of (61) has been given, to sixth order, by 
(46) or approximately by (48) «. The second term, the theta dependent term, is given 
approximately by: 


^ (0) - g" sin i sin^2 0 + ^ sin^i sin^2 0 


(62) 


n 


n 


To fourth order the total sensitivity is then: 

2 4 

5X „ q sin i 1 sin i ,, , . 2„ ^ 

X 4 + 8 “4““ (1+c l sin 20) • 


(63) 


At this point some insight into the angle sensitivity characteristics of the 
WFB element can be obtained. Clearly the second order expansion, equation (19), 
is no longer valid when the second and fourth order terms are equal in magnitude. 
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This occurs when. 


2 

sin i = 2 | q j n 


i = sin 1 (2 |q|n^) 


1/2 


C64) 

(65) 


For this angle the wavelength shift, to second order, would be 



2 

q , 


( 66 ) 


Since filters or spectrometers are normally used with the maximum wavelength shift 
equal to the full width at half maximum (FWHM) , the second order expansion should 


be valid for filters with FHWM less than the limit given by equation (66) . 

For quartz and calcite, the most commonly used birefringent materials , i and <5 

e 

at 5000& are (9.7 > .08799^) and (45 , 34.7^) respectively. The second order 
approximation is, therefore, almost always valid for wide field calcite elements, 
but is often not valid for quartz elements « In fact , when a quartz element is used 
beyond 9.7 even the sign of the wavelength shift predicted by the second order 
approximation is incorrect. 


X 


It should be noted that all the terms in the expansion for A (equation 42 
et. al.), with the possible exception of the second order term, are positive, 
Since the coefficients do not decrease rapidly, the overall convergence of the 
series is, therefore, slow. Detailed calculations of wide field systems must 
consider terms of higher order than the fourth, especially when the second order 
term is negative. 
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V. MISALIGNMENT AMD MISMATCH EFFECTS 

The variation of the angular sensitivity with azimuth, theta, is negligible 

for well aligned WEB elements. If the crystals are not perfectly orthogonal, 

however, then the situation changes markedly and the field of view develops a 

* 

strong azimuthal dependence. It is worthwhile to investigate these effects of 
non-orthogonal crystals not only to determine the error sensitivities, but also 
because angular misalignment offers a mechanism for cancellation of the second 
and fourth order terms, thereby extending the field of view for some angles. 

When the two crystals in the WFB element are non-orthogonal the time delay 
is again given by equations (.39) through (43) but with 

2 n 

■z f = sin (0+90+a) = cos (0fa) £67) 

in which a is the misalignment angle* The time delay is again expanded as a power 
2 

series in x , the angle of incidence dependency , but the polynomial factors, 0 (n) , 

WF 

need to be modified. In place of equation (51) the expression for 0 (n) now 

WF 

becomes: 

n 

= 1 + h C-l ) 1 (1-pV (?) Csin 2i 0+ COS 21 (0+ o))n>l. (68) 

i » 1 1 


By analogy with equation (52) we will define 

(a) = sin 2 *0 + cos 2 * (0+ a). 

Letting 

2 1 

R « sin a cos 2 0 + — sin 2a sin2 0 > 
we have 

2 2 

cos (-0 + a) ~ cos 0 - R 


(69) 

(70) 


. (71)a. 
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and 


• x 

cos 21 ( 0 + a) = cos 2 l 0 + x (~l) k (£) R k cos 2 ^“ 2k 


(71) 


The factors, S ± (a), may then be expressed as the well aligned component, S ± , 
plus a correction sum: 

S i <“> - s i + E <-D k O E k .<=os 21 - 21 & . (72) 

k = 1 

Thxs means that the polynomial factors, O^W, may also be decomposed into a 

well alxgned component plus an alignment error component- The alignment error 
component has the form: 

” 2 1 £. ! C_1) (1 ' P 5 ( i> k L (- 1 ? <P R k cos 2i ^ k 


(73) 


For small misalignments, the factor R may be approximated by 

R - a sin 2 9 , ( 74 ) 

and terms in R or higher powers of R may be neglected. Equation (73) is then 
replaced by: 

Cn;a) _ a s in 2 0 £ . -,i „ 2 ,i , n% 


0 . 


WF 


2 q 


L (-D 1 ( 1 -P 2 ) 1 (?) i, 
i = 1 1 


(75) 


which equals 

na sin 2 0 p 2n_2 

Again letting p be very near unity. 


(76) 


^WF^ n,a ) ~ sxn 2 0 v 


(77) 
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Adding this terra into equations (57) yields, respectively: 

0^(1) = ^ + a sin 2 0 

0 W pC2) = - 1 - q sin 2 2 0 + 2a sin 2 0 

, N 2 (78) 

Oyp (3) - - 2 - 3q sin 20 + 3a sin 20 

2 3 t 

°WF^ - - 3 - 6q sin 20 + q sin 20 + 4a sin 20 . 

Including the misalignment terms, the angle sensitivity relation is to fourth 
order : 


6X 

X 


, - 2 * 

±. sm i 

4 2 

n 

o 


. 4 . 


(q + 2a sin 20) + j 1 (l-2a sin 20+ q s in 2 20). (79) 

n 


Developed this -way it is easy to see the relative magnitudes of the theta independent 
component, the well aligned theta dependent component, and the misalignment theta 
dependent component. As may be seen, an angular misalignment can be chosen to cancel 
the other components for certain directions. This will be discussed further subse- 
quently, but first we will show a similar effect for a length mismatch between the 
two halves of the WFB element. 


Suppose that the lengths of the crystal parts are 
d 




+ Td, 


(80) 


d 2 2 Td 'o 


The time delay is then 
& 

A = ^ 


+ TA 


f (x, p 


, q, z) + f (x, p, q, z ') I 

1 J 


f (x, p, q, z )- f (x, p, q, z ’) 


(81) 


The first term is the familiar wide field time delay, while the second term is 
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similar except for the minus sign. Again, this second term may be expanded as 

2 ' 


a power series in x ,but with- the polynomial factors: 

°WF (n;T) = 2q [ (1 - (1 ~P 2 > z) n - a-Cl-p 2 ) z r ) n J . 

Xf we consider only a length, mismatch, this becomes 

n 

°WF^ n ’ T ^ ~2q (“1) (l-p 2 )* (J) Csin 2i 0 - cos 2i 0-). 


C82) 


(83) 


The factors 


Ih = sin^G- cos^G 


C84) 


may be expressed as 


i-1 

- 2 * 


D = - cos 20 


f- £ 


k - 1 


“I 

B . (k) sin 2 G 
1 . 


(85) 


in which the coefficients, (k) , may be determined from the relations: 


B. (1) - ifS 


i>2 


( 86 ) 


B i (k > = "4 B i-2 Ck_1) ~ A^iCk) i>3, k>l 

The terms, A.^_^(k), are a s defined by-equations (54). 

Since we are assuming that the length mismatch, T, is small and that p is 
near unity, we may neglect all but the first term of (83): 

0 WF (n;T) *s n cos 2 0 . (gy) 

The contribution of the length mismatch to the time delay expansion is, therefore. 


A (T) A q T cos2 0{c‘^x + 2c 2 x. +. . ,+nc^ x 2n +...} (88) 

Combining this with the fourth order expansion, equation (63) , thje total sensitivity, 
including mismatch, is: 

2 4 

■ ii - - 4 SX 2 1 <9 + 21 ct>s 2 0) + ir — ” 1 u + q sln 2 2e - 21 cos 20 ) . (89) 

A » V 8 -o 4 
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Equation (89) shows that the effect of a length, mismatch is very similar to a 
rotation misalignment except that the field of view expands in different angular 
ranges . 

We will not consider the combined effect of rotation misalignment and length 
mismatch since the additional terms would be on the order of (aT) , which can be safely ■ 
neglected. 
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VI . THROUGHPUT ADVANTAGES OF WIDE FIELD SYSTEMS 

With the relations between relative wavelength shift and the angle of 
incidence and azimuth derived in the sections above it is possible to demonstrate 
the relation of the WFB element to the F— P, MI, and WFM. For negative crystals 
there are two regimes as indicated by equations (.64) through (66) , depending upon 
whether the second or fourth order terms dominate the angle sensitivity. 

When the second order term dominates 

ii ^ q 2 / 2 , (90) 

X 

and the maximum incident angle for an allowed relative wavelength tolerance is 
obtained from 


’ 2n 

sin i = o 


— A' 

1/2 <■ X ' 


Relation (9) is a similar result for the F-P and MI. So that the ratio of acceptance 


angles of a WFB to a F-P or MI is: 


* ■ £ 


1/2 . 


Since q is always less than one, the WFB has a larger field than the F-P and MI. 

The gain in throughput for equal area devices, neglecting vignetting is 
2 ? 2n 9 

„ T “ 2. / O N „ V 


g 2 - V 


( — o y 


For a calcite filter (q - .12) the gain in throughput is approximately 46 for 


n Q ** 1.67 and n pp = 1. 


Regardless of the sign of‘ q, when the fourth order term dominates the wavelength 
sensitivity the maximum incident angle for the WFB is obtained from 
sin i - 2 n o <| ^) lM = n Q (8 ^) lM 


( 94 ) 
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and the ratio of acceptance angles with- respect to a F-P or M is 


i, - (r*> a 


1/4 


n 


FP 


SX ; 


(95) 


and the gain in throughput is 
n 2 1/2 

FP (96) 

<5.X 

Since — can be very small the WFB has the potential of significant throughput 
gain with respect to the F-P and MI. Comparison of relations (94) and (25) shows 
that the angular field and hence throughput of the WFB and WFM are equivalent when 
the fourth, order terra dominates. , 

Por positive crystals another possibility exists because the wavelength 
sensitivity is zero at the point the second and fourth order terms balance in 
magnitude, thus the wavelength sensitivity has an extremal which is determined by 


3_ ( Qs _ 0 

9 (sin i). ~ ° 


(97) 


and occurs when 


sin i = n Q ^ q . 


(98) 


From relations (98) and (63) the extremal value is 


~ - - q /S . 


(99) 


From the equality condition (equation 64) the relation between angular field and 
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wayelength. tolerance in the optimized condition is 


sin t = 2 n (2 B-) 1 / 4 
o J 


( 100 ) 


So that the gains in field and throughput compared to a F-P or MI are 


n o 1 X 1/4 

I - 2 fcp ) (f i ) , 

° n FP ^ 


( 101 ) 


and 


n 0 2 1 X 

G 0 * 4 (-^ ) 2 (f | ) 3 

° n FP 2 oX 


( 102 ) 


respectively* The optimized condition, thus yields a factor of 2 gain in through- 
put over the fourth order condition. It should also be clear that a WFM can be 
adjusted to increase its throughput by an appropriate adjustment. 

The angular characteristics of the WFB can be illustrated by a pair of log-log 
plots as shown in Figure 3. The upper plot, (3a ) 9 shows sine i versus relative 
wavelength, dA/X^, for quartz and calcite. The true values (solid lines) are 
compared with the 2nd and 4th order approximations (dashed lines). In addition, 
a curve for a crystal of optimal q (equation 100) is also presented for comparison. 

On the right hand scale of Figure 3a are values of the angle of' incidence and on 
the top scale are the values of filter bandpass at A q = 5000& for telecentric systems 
whose maximum shift equals the wavelength tolerance. Immediately below, in Figure 3b, 
are the corresponding plots of the ratio of incidence angle (left scale) and the 
throughput ratio (right scale) versus the relative wavelength tolerance. 
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Figure 3 clearly indicates the regions for which the 2nd and 4th. order ap— 

proxitnat ions are valid for quartz and calcite* Xt also presents the size of the 

gains in acceptance angle and throughput which are possible with a WFB with respect 

to a F“P or MI, A calcite WFB can be built for the range of relative wavelength 
-5 -3 

tolerance from 10 to 10 where a throughput gain of nearly 50 can be realized. 

For quartz, the wavelength tolerance range from 2 x lO - ^ to 6 x ICT 3 is practical, 
for which throughput gains from 40 to 300 are possible. 

The arguments above for optimimization of q are not strictly academic. Although 
the selection of natural crystals is strictly limited it should be possible to create 


sheet plastic materials with, desired q. Uniaxial sheet plastic wide field elements 
can be built by lamination. 

Even if it is not possible to select q it is clear from equations (79) and (89) 
that the factor multiplying the second order term can be made zero by appropriate 
rotation or mismatching of crystal lengths, respectively. 

Adjustment of the angle or lengths has the greatest effect when the second order 
term dominates the wavelength sensitivity. To estimate the field gain consider 
equation C79). and suppose a is adjusted such, that at 45° and 225° the second order 
term vanishes. Then the wavelength sensitivity function will have an oval shape 
with major and -minor axes given by 


sin i 

major 


2 n & 


1 6X 


1/4 




C103) 


and 


sin i 


minor 


-#-n o 


(104) 
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Th.e misadjustment results in an increase in throughout 
I | 1/2 1/4 

w - ¥ aos) 

compared to the properly aligned element. For a calcite element of O.lS bandpass 
(X500 o£) this gain is about 3, which means that the field in one direction in- 
creases by a factor of 6. "Using the higher order terms yields somewhat higher 
estimates. 


A-34 



- 28 - 


VII. FRINGE PATTERNS OF WIDE FIELD SYSTEMS 

An important use of the wavelength, sensitivity calculations performed 
above is for the explanation of the fringe patterns produced by wide field ele- 
ments of different materials in various states of alignment or length 
mismatch. The patterns are easily produced by observing an element over a diffuse 
monochromatic source. Since a bright fringe occurs whenever the optical path 
dxfference is an eyen multiple of halfwaves, the fringe pattern is a contour map 
of the wavelength, yariat ion. That is, the radius vector to a given bright fringe 
is determined by the angle of incidence and azimuth required for a particular wave- 
length, shift. It also should be noted that unless the illuminating source has two 
wayelengths the contours are unsigned; the fringes can be caused by either a posi- 
tiye or negative wavelength shift. Of course, with two wavelengths reversal of the 
wavelength shaft wall cause a reversal of the fringe colors 

Once the. fringe patterns are understood "the observer has a wonderfully sensi— 
tpye tool to use for aligning and matching crystals. Further, the patterns are very 
beautiful and are in our opinion worth understanding on purely esthetic grounds. 

The bright fringes occur from the fourth, order expansion whenever 


- 2 - . 4. 

m = R f - Sln - 1 1 Cq + 2a sin 20 + 2T cos 20) + i 

4 


L 4n 


8n 


Q06) 


where m is a signed integer, and 



o 

For equation (106) it has been assumed that 
4th order term in sine i can be neglected. 


(107) 

the effects of q, a, and T on the 
In order to calculate the fringe pattern 
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accurately, higher order terms must be considered in regions where the wavelength 

sensitivity is slowly yarying. Contours on gently slopes can move a great deal 

with very small changes in inclination of the slope. For this reason a computer 

program called PAT has been written to solve an eighth degree polynomial equation. 

Since all terms higher than second order in sine i are positive (equation (46) , 

the general shape of the fringe pattern can usually be determined from equation 

(106), which means most of the properties of the fringe patterns follow from the 

2 

characteristics of a simple equation in x 


4 

.x 



C108) 


where 

q* = q + 2a sin 2 0 + 2T cos 20 . 


(109) 


Fringes arise from solutions of relation (108) for all possible integer values 

(both positive and negative) of m for fixed q, a, and T. The solutions are limited, 

2 -2 

however, by the condition that x be real, positive and less than n . The limi- 

o 

2 

tations on x are because sine i must be less than one. As is well known, there 
are two solutions of a quadratic equation' which involve, respectively, a positive 
and negative contribution from a surd. For equation (108) the surd is 


Q “ (2q f ) 2 + 32m/R 

- 


1/2 


Because x must be real Q must also be real. 
2 

that is, when q* < 8m/R • 


( 110 ) 

This means that when q 1 is small, 

( 111 ) 


only solutions for positive m are possible, while if q' is large 

8m 


i . e . , q ’ " > 


R 
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both positive and negative solutions are possible. The condition that x be 
positive places further limits on m depending on whether q is positive or negative. 
The valid solutions for positive and negative in with the assumption that either 
(111) or (112) holds are summarized in Table I. 

For perfect elements (a and T equal zero) of small q, Table I shows that 
there are only solutions for positive m regardless of the sign of q. Thus the 
fringe pattern for small q is generated by 


sin i = n 

o 


,8m. 1/4 
t R )■ 


(113) 


Note that the circular fringes' have radii proportional to the fourth root of the 
integers and, therefore, form a pattern similar to that of a WFM. 

For perfect elements of large negative relative birefringence, Table I shows 
that there exist only solutions for positive m and that the fringe pattern is gener- 
ated by 


sin i = 2n 

o 



(114) 


The fringe radii of this pattern are proportional to the square roots of the inte- 
gers, and thus the pattern is similar to that of a F-P or "MI. 

For the case of large positive values of q. Table I suggests that there are 
three possible sets of fringe patterns. For negative m there is the expected set 
of fringes given by 


sin i = 2n 

o 


(M ) 1/2 

^Rq } , 


(115) 
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and also the set given By 


sin i - n 

o 


2q (1 + 



(116) 


Further, there are positive ra fringes giyen by 

b * - ^ 1 m 


sm 1 = n 


(117) 


The fringes indicated by expressions (116) and (117) occur because of the cancel- 
lation of the second and fourth order terms, and the angle 

i = sin -1 (2qn 2 ) 1 / 2 . 

Q (118) 


represents the transition between second and fourth order fringe patterns. 

Unfortunately, positive crystals with high relative birefringence are not 

normally made into wide field elements, so that these fringes have not been observed, 

at least by our group. However, it is straightforward to calculate the appearance 

of any fringe pattern using PAT. Shown in Figures 4a, 5a, and 6a are the locations 

of the first few bright fringes of a (A 5000, ^ ~ 2 x 10 wide field element 

6 

of quartz (small q) , calcite (large negative q) and reverse calcite, (large positive 
q) , respectively. In all PAT plots negative fringes are dashed. Reverse calcite is 
imaginary material with reversed ordinary and extraordinary indices from calcite. 

For elements that are misaligned, a is non zero and the value q 1 oscillates in 
magnitude when 2a is less than |qj and oscillates both in magnitude and sign when 
2a is greater than |qj. First consider the case of a large negative q material. 

When 2a is less than |q|> |q T | will decrease in the 45°-225° direction and increase 
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Figure 5. Fringe patterns for misalignments of 0, .2q, .4q, .5q, . 75q, 
respectively, of a calcite element. 
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in. the 135 —315 direction. Since the fringe radii are inversely proportional to 
[q [ , the fringes expand along the 45 —225 direction and contract along the or- 

thogonal direction. The expansion and contraction are, therefore, along and or- 
thogonal to the polarizer direction. When a is positive the expansion is along the 
45°-225° direction, while when a is negative it "is along the 135°-315° line. Note 
also that as a approaches jq| the fringes do not extend to infinity, but are limited 
by the small q solution which is almost independent of q and depends mostly on m and 
R. 

When 2a is greater than |q|. Table I, shows that fringes of negative m become 
possible. The negative m solutions give rise to egg-shaped fringes that nest inside 
the extended positive fringe system. The fringe pattern for a l£ calcite (large 
negative q) WFB element is illustrated in Figures 5b, c, d, e where a equals 
•2q, . 4q, ,5q, .75q, respectively. The pattern of Figure 5e is very similar to that 
of a biaxial crystal. 

For misaligned WFB elements of large positive q the situation is somewhat 
different. When 2a is less than q the inner fringes behave in a similar manner to 
those of a negative q material except the elongation is along the 135°-315° line 
and the contraction along the 45°-215° line for a positive rotation. Both the 
positive and negative m outer fringes expand along 45°-225° direction and contract 
along in the 135 —315 direction because (see Table I) they are proportional to q * 
rather than inversely proportional to j q 1 |. as are the inner fringes . When 2a is 
greater than q the angular range of the negative m fringe system is limited because 
there are no solutions for both negative to and q'. What occurs is illustrated in. 
Figures 6b, c, d, e for a 1& reverse calcite (large positive q) WFB element. The 
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a values are the same as for Figure 5. 

When q is small the fourth order term in sine i dominates and quite a differ- 
ent process dominates the rotational misalignment fringe patterns. This can be 
seen using equation (79) and the approximation that q can be neglected with respect 
to a. We then have: 

4 ? 

m/R - (1 - 2a sin 2'0) x - 2a sin 20 x , (119) 

When a is small compared to 1/R only positive m solutions are possible and the 
fringes will elongate along the line 45°-225° and contract along 135°“315°* If 
a is increased further a negative fringe system .will occur inside the elongated 
positive fringes . The result for small q is that regardless of sign the fringes 
will elongate along the 45°-225° or 135°-315° line depending upon whether a is 
positive or negative* The first negative fringe -can be expected when a is ap- 
proximately equal to: 

a Q = 1 /R . (120) 

The value a Q in (120) is determined by the extremal value of equation (110) 
for the case m equals minus one. Shown in Figures 4b; c, d, e are the fringe 
patterns of a quartz 1$ WFB element for a equal .4 a Q , .8 a Q , a Q , and 1.5 a Q . 

The patterns of Figures 4, 5, and 6 represent a reasonably complete catalog 

of the misalignment fringe patterns of uniaxial WFB elements. When WFB elements 

are mismatched in length a very similar set of patterns arise as should be expected 

'from comparing equations (79) and (89). The elongations are along the 0°-180° or 
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90 -270 lines or along or orthogonal to the optic axes depending upon whether T 
is positive or negative . 

In general a WEB element if it is remoyed from a birefringent filter or in 
an early state of manufacture will have both angular and length mis adjustment ♦ 
Therefore, the field of view may be elongated in an arbitrary direction. By care- 
ful adjustment of the element halves with respect to each other, however, it should 
be possible to get either circular fringes or more likely fringes elongated along 
or orthogonal to the optic axis. Erom Eigures 4, 5 or 6 and equation (89), the 
elongations can be estimated. Nearly perfect fringes will occur for large (small) 
q WFB elements when a and T are less than 1/20 q (1/10 R -1 ^ 2 ) . For calcite and 
quartz this means the tolerance on angle and length match of a l£ are (.32°, .004 cm) 
and C-114 , .026 cm), respectively. 
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VIII. DISCUSSION 

In this paper the main intent has been to explain the very significant 
throughput advantages of WFB filters and the possible modifications of the 
field of view for sources that are extended in one direction. The paper is 
not intended to be a complete discussion of off axis error sources and light 
loss in wide field elements as was paper III 8 of this series for on-axis effects. 
The effects of imperfect waveplates, waveplate field, and crystals with optic 
axis not in the plane of the crystal face will be considered in a further paper 
of this series. * 

This paper has resulted in large part from experience gained in constructing 

three wide field of view Birefringent filters. The filters are a pure Lyot 10? 

bandpass quartz filter, a redundant alternate partial polarizer (RAPP) 60& bandpass 

polyvinyl alcohol filter, and a modified alternate partial polarizer filter (APP) 

with 60m£ bandpass. All of these filters use the plastic waveplates discussed in 

9 10 

paper I through III ’ and the internal apodization of pulse structure concepts 
introduced in paper IV . However, a new general technique for partial polarizer 
filter design was used for the RAPP and modified APP designs which also will be 
described - in a later paper of this series. 

The quartz and polyvinyl alcohol filters are wide field and tunable from 
4500-7000A and have diameters of 70m and 35mm, respectively. They exhibit the 
large throughput advantages predicted above. The polyvinyl filter is especially 
interesting because its elements are constructed by laminations of one wave, sheets 
of stretched polyvinyl purchased from Polaroid Corporation. Such lamination tech- 
niques should allow experimentation with optimally stretched sheet plastic materials. 
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GENERAL INTRODUCTION 


This reporx is on work done for the NASA Solar Physics Spacela’o 
Quick Reaction and Special Purpose Facility Definition Team on the critical pro- 
blems of a space qualified birefringent filter suitable for use on an Ha/magnetograph. 
This report uses data , fabrication techniques and computer programs for the analysis 
of birefringent systems previously developed by the Lockheed Solar Observatory. 

The construction design philosophy of the filter modules tested as well as the 
construction techniques were based on previous experience with numerous birefringent 
filters. 


The critical problem for flight of a birefringent filter is the shock 
mounting of the calcite. The design presented here bonds the calcite block 
with silicon rubbers to the calcite holder. The calcite together with its 
all necessary polarizers and rotating achromatic plates are mounted together 
m units called a 1 filter module. By using a set of modules containing calcite 
crystals of differing lengths, a filter can be produced. A description of the 
modules is contained in Section I. Also described in Section I is a container 
for the filter modules, which can be used both to hermetically seal the system 
or contain an index matching oil. 

A second serious problem with birefringent filters has been wavelength 
drifts introduced by thermal variations. The traditional approach to the thermal 
drift problem has been a hardwired temperature controller and oven enclosure 
for the calcite. Section II reports on the response of a filter element while 
being controlled by the Lockheed Temperature Control. 

A tunable filter such as proposed for the Ha/magnetograph offers the 
possibility of using the tuning elements to compensate for any thermal drifts 
thus eliminating any need for a hardware temperature controller if the rate of 
change of wavelength with temperature is known. Section III describes the 
determination of the wavelength sensitivity to temperature of calcite. The 



Lockheed. Universal filter has been built with temperature sensors on The 
calcite elements. Section XV describes operation of the filter using a software 
control algorithm instead of a hardware Temperature controller. 


SecTion V briefly reviews some radiation considerations of : 
terns. Also included is an appendix That indicates future areas of 
for filter maTerials and systems. 


filter sys- 
researeh 
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I. SPACE HARDENED FILTER MODULES 


INTRODUCTION 

The goal of the space hardening experiments was to contain calcite elements 
in a cell that would sustain both shock testing and thermal cycling. It was 
desired that the cell also contain the necessary components for a complete 
tunable module. A preliminary test model was used to verify the cell design, 
a prototype was designed to verify the tunable module, and an oil sealed 
prototype to verify a complete filter assembly. 

In addition, optical cementing techniques were developed to allow reliable 
cement bonds between calcite to calcite, calcite to glass, and calcite to 
polyvinyl alcohol (PVA) . Bonding was of some concern because previous 
calcite cements had a tendency to cause fracture when bonded assemblies were 
thermally cycled. Reliable cements allow construction of filters that do not 
require index matching oils. 

SHAKE TEST 

Preliminary Test Unit 

Before making a f u ll 30 mm aperture 1A filter module, a test unit; with 
25 mm aperture was constructed. The test unit consisted of 4A and 2A wide 
field elements. Table I lists the components of the Mark I test unit. The 
test unit was constructed with and sealed in its housing with GE RTV 602. 

The test units were shaken to the specification given in A.stronomy Spacelab 
Payload A.ccommodations and Interfaces Handbook (May 1975), Figure 4. 2. 2. 2-1 
(reproduced in Figure l). The unit was shaken in three orthogonal axis, 
where the y axis was parallel to the filter optical axis. A sample shake 
test chart is shown in Figure 2. 

The first rest unir shake, 2-25 - 7^, showed separation of the waveplates 
and rheir cover glasses. The separation was evident from fringe patterns 


B 3 



seen when the filter was mounted over a diffuse helium source. A photograph 

of the -fringes is shown in Figure 3. The filter was disassembled and the 

separation of the waveplates was verified. The failure was attributed to the 

poor bond between RTV 602 and the polyvinyl waveplate material. 

♦ 

Although the separation occurred, the calcite showed no evidence of 
damage nor did any of the components rotate with respect to each other. 

The spectral performance of the Mark I version of the test; unit did not; 
change. 

The test unit was rebuilt using Eastman Kodak HEF-4 optical cement -co 
bond the waveplate material. The glass covers of the polyvinyl in intermediate 
positions in the filter were eliminated. The optical component lisp is 
shown in Table II. The filter was bonded in the cell with RTF 602 as before. 
The mechanical drawing of the Mark II system is shown in Figure 4. A photo- 
graph of the cell and shake fixture are shown in Figure 5 . 

Kodak HEF-4 is a non-hardening optical cemen-c . Calcite elements have 
been cemented with HEF-4 and temperature cycled between 20° and 60°C. 

On 3-17-76 the unit was reshaken. No separation was observed. No 
evidence of any adverse effects was shown by the filter. Spectral performance 
was not affected. 

Prototype Filter Module 

After the successful experience with. the test unit a flight prototyue 
unit was constructed. The flight prototype unit has an achromatic half wave- 
plate and a rotating polarizer. Its optical components are listed in Table 
IH. A mechanical diagram and a photograph are shown in Figures 6 and 7 . 

On 5-4-75 the prototype unit was successfully shaken. The calcite did 
not rotate and the spectral performance was unchanged. 

The calcite used was fabricated by Halle Nacht and has two small cracks 
on its edge. This calcite was especially chosen to see if the cracks would 
increase. They did not. 


B- k 



TABLE I 


MARK I COMPONENTS LIST 


Material 

Thickness ( cm) 

Orientation 

Comment 

1. 

Glass 

.203 


Bonded -do 2 

2. 

1/4 X (FVA) 

.00254 

^5° 

3 

3. 

Glass 

.203 


lx 

4. 

Calcite 

.254 

0° 

5 

5- 

Glass 

.203 


6 

6. 

1/2 X (FVA) 

.00254 

VJ1 

O 

7 

7- 

Glass 

.203 


3 

8. 

Calcite 

.254 

90° 

9 

9. 

Glass 

.203 


10 

• 

O 

H 

EN-38 (FVA) 

.00254 

0° 

11 

n. 

Glass 

.203 


12 

12. 

Calcite 

.127 

90° 

13 

13. 

Glass 

.203 


14 

14. 

X/2 (FVA) 

.00254 

45 0 

15 

15. 

Glass 

.203 


16 

16. 

Calcite 

.127 

0° 

17 

17. 

Glass 

.203 


18 

18. 

X/4 (FVA) 

.00254 

45° 

19 

19. 

Glass 

.203 
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TABLE II 



MARK II COMPONENTS LIST 


Material 

Thickness (era) 

Orientation 

Comment 

1. Glass 

.203 


Bonded to 2 

2. 1/4 X (FVA) 

.00254 

45° 

3 

3. Calcite 

.254 

0° 

4 

4. 1/2 X (PVA) 

.00254 

45° 

5 

5. Calcite 

.254 

90° 

6 

6. HU-38 (FVA) 

.254 

0° 

7 

7. Calcite 

.127 

O 

O 

o\ 

8 

8. X/2 (PVA) 

.00254 

45° 

Q 

y 

9 . Calcite 

.127 

0° 

10 

10. X/4 (FVA) 

.00254 

^5° 

11 

11. Glass 

.203 




EEF4 
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Material 


TABLE III 

Prototype Components List 


Comment 


Thickness (cm) Orientation 


1. Glass 

.203 

Rotates 

Bonded to 2 with HEP -4 

2. M-38 

.00254 

as 

3 

3. Glass 

.203 

unit 


4. Space 

.114 



5. Glass 

.203 


Bonded "&o 7 with HEF-4 

6 . X/4 (FfA) 

.00254 

45° 

8 

7. Calcite 

.5588 

0 ° 

9 

8 . X/2 (PVA) 

.00254 

15 6 

10 

9. X/2 (PVA) 

.00254 

75° 

11 

10. X/2 (PVA) 

.00254 

15° 

12 

11 . Calcite 

.5588 

0 

O 

On 

13 

12. M-38 

.00254 

0 ° 

14 

13. Glass 

.203 
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SEALED FILTER ENCLOSURE 


After the successful shake test of the prototype filter module the spur 
gears were replaced by bevel gears which were more suitable for the rotary 
feed-through system used for the sealed enclosure assembly. The enclosure 
was designed as a hermetic or oil sealed case for the 'oirefringent filter. 

The enclosure actually built is just sufficiently long to contain the prototype 
filter module, but a larger enclosure with more feed-throughs could contain 
any number of filter modules. 

The enclosure was built with 0- ring seals at all the non- rotating joints. 

The enclosure components and the complete enclosure are shown in Figs. 8 and 9. 
For a flight unit epoxy or welded seals might be judged more suitable. The 
rotary feed-through necessary to control the tuning module has no rotary oil 
seal. Rather, rotary motion is transferred through an intermediate plate 
that rocks. The plate is welded to a metal bellows which in turn is welded 
to the side of the feed-through so that there is no piercing connection between 
the input and output shafts. No magnets of any kind were used in the seal. 

The enclosure was tested for gas tightness by flooding the exterior with 
helium while pumping on the interior with a vacuum pump. A helium leak 
detector in the output part of the pump showed no leakage signal. The tuning 
motor was operated as the unit was pumped. Following the gas seal tests 

the enclosure was filled with index matching oil and placed in a vacuum chamber. 
After lorty hours in the chamber at 10 °torr no oil was detectable on the 
surfaces of the container. The motor was operated periodically while the 
unit was in the vacuum chamber. 

On 6-25-76 the oil filled enclosure was shaken to the shuttle shake profile. 

At the time of the shake ~ est the enclosure contained the prototype filter 
moaule. No oil leakage was observed after the shake. The unit was then reulaced 
in the vacuum chamber and left in vacuum for 16 hours. The motor was periodically 
rotated. After removal from the vacuum chamber, no l eakag e was observed. 



The oil inspection was carried, out by wiping a clean cotton swab over all 
metal parts. The exposed optical surfaces were visually inspected. Any 
oil filter would create interference colors or tint, none were observed. 

A-R COATING- POSSIBILITIES 

A marjor problem with ground based filters has been leaking and contamination 
of the index matching fluid or oil in the filter. In space, the consequences o 
leakage are much more severe. One solution to the leakage problem is to 
completely seal the filter as described above. Another is to eliminate the 
index matching fluid. 

The index matching fluid serves two purposes . The first is to eliminate 
fresnel reflection losses at the interfaces. The second is to eliminate 
image displacement caused by rotating tilted optical elements in the filter. 
During this study a cementing technique has been developed which allows 
reliable stress-free assembly of calcite to calcite and calcite to glass. 

The index of the cement is 1.57. The index of the cement used for the 
waveplate is 1.54. 

If a Lyot filter is constructed from a set of turnable modules of the form 
half wave, quarter wave, calcite, half wave, calcite, polarizer where the 
first half wave plate rotates and all the other components are cemented 
together, there are only four air-glass interfaces per module. A plot of 
maximum transmission, neglecting polarizer losses, versus number of modules 
is shown in Pig. 10 for several values of per surface reflection loss. 

Shown in Fig. 11 is the per surface reflection loss for a s amp le four 
layer antireflection coating. Prom 4300A to 6600A the per surface 
reflection loss is less than 0.2$. This spectral region spans the useful 
region of currently available polarizing materials. Recent developments 
suggest coatings may be available that have reflectance less than 0.2*$ 
over the entire range from .3 to 2.2 microns. 

For the filter design suggested above, only the half wave plates rotate. If 
the half wave plates are 0.25 cm thick and are mounted in 5-cm ring 

_3 

bearings with less than 1.25 x 10 cm axial play, then the displacement of 
the optical beam through the plate is 
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d 


d = (2.5 x io" 1 ) ( 1,2 \y 0 3 ) (xt||) 

d = 4.14 x 10" 5 cm 


An image scale of 40 arc seconds per cm is such that 0.1 arc second occupies 
2.5 x 10 3 cm, the size of a typical diode array element, a *j-0 tt cm scale 
thus is nearly 2 orders of magnitude greater than any image motion caused 
by rotation. 

Tb verify that a cemented filter would operate well, one of the observatory 
Halle H a filters was disassembled, degreased, and reassembled with cemented 
components and A-R windows . The filter is now in routine use in the 
observatory patrol system. 

The rebuilt Halle has been repeatedly cycled between 10° C and 50° C. At 
present it is left outside with the power off during the evening hours and 
is also uncovered during normal operation in order to subject the cement 
bonds to larger than normal thermal loads. 


B- 10 



TEMPERATURE COIWROIiLER TESTS 


PURPOSE 

The purpose of the thermal controller tests was to evaluate the ability 
of a stand alone temperature controller to maintain the temperature of cal- 
cite filter elements in the presence of environmental thermal loads and 
solar heating. The temperature controller used was the Lockheed Filter 
Controller designed by R. Lindgren. This device senses ambient temperature 
and supplies a steady state heat input based on the difference between ambient 
and the control temperature. Xn addition, it applies a fine control heat 
input based on a sensor near the calcite crystal. The power input follows 
the law: 

HI = k 1 (t c -t a ) + k 2 (Vt s ) 2 sgn(t c -t g ) + k 3 

where k^, k^ and k^ are adjustable parameters, ana t^., t^ and t are respective 
the control, sense, and ambient temperature values. If the sensed temperature 
is greater than the desired control temperature, the ambient control heat 
input is reduced. The control input is limited to be less than 20 percent 
of the ambient heat input. 

M&TH0D OF MEASUREMENT 

Two fundamentally different systems were used to measure the characteris- 
tics of the controller. One, using a quartz crystal thermometer, measured the 
temperature in the oven close to the calcite. The second measured the mean 
temperature through the crystal using a laser probe. The laser probe allowed 
measurements of possible temperature inhomogenities in the calcite it- 
self. Measurements were made on a .081 wide field filter element which was 
installed in a Halle oven. A quartz thermometer probe was mounted in the 
aluminum block enclosing the calcite (see Figure 12). The temperature was 
read out on a Hewlett-Packard (Dymec) Quartz Thermometer. The quartz 
thermometer is repeatable over the short term to ± ,0005°C and over times 
of days to ± .001°C. It is linear over the range used to ± .0008 degrees 
per degree. 
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Shown in Figure 13 is an optical schematic of the test system. A 
2 mm He-Ne collimated laser beam is sent through the filter and xhe output 
detected by a photometer. In addition to the laser beam, a solar beam can 
be simultaneously shone through xhe filter. The filter is tunable by rotating 
a half waveplaoe behind the entrance polarizer. Both the entrance and exit 
polarizers are mounted separately from the Halle oven. 

As the temperature changes, the index of refraction difference and 
the leng-h of the calcite changes, so that the laser line is alternately 
transmitted and extinguished. For this crystal there are approxim at ely 2.22 
transmission cycles/°C. Fractions of a cycle can be measured by rotating 
the half waveplate to an angle where a transmission minima occurs. The 
locations of transmission minima can be reliably located to one degree in 
angle, thus the equivalent temperature measurement accuracy is ± l/90 of 
1/2.22 or ± .005°C. 

CHARACTERISTICS OF THE CONTROLLER 

Four identical controllers were evaluated. Controllers were set to con- 
trol at approximately 30°C, 36°C, 42°C and 48°C. In a normal room environ- 
ment, all the controllers would maintain an oven temperature of ± ,001°C 
as indicated by the quartz thermometer readout. When the temp'erature reading 
of the quartz thermometer was stable, the laser transmission was also stable. 

The effects of thermal environment of the temperature controller its elf 
was tested by separately placing the calcite oven and the temperature control 
electronics in a thermal chamber that could be cycled from 0° to 50°C. 

The temperature of the calcite oven changed by - . 0028°C/°(J of ambient over the 
ambient range from 0°C to 50°C. The electronics caused a change of .001°C/°C 
ambient. With precision components the temperature control system could be 
improved by a factor of at least two. 
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SOLAR EFFECTS ON OVEN TEMPERATURE 


With, the oven interior temperature well stabilized, ana the waveplate adjusted 
so "hat a transmission minima occurred for the laser beam, an f/l? solar beam was 
allowed to pass through the filter. The beam was prefiltered by a 8 a three oeriod 
blocker and an IR rejection filter. There was some tendency for the trans- 
mission minima to move a degree in angle or for the calcite to change by .005°C 
in temperature. However, after repeating the test a number of times, no 
statistical significance could be obtained for a shift of a degree. 

To further test for evidence of solar heating, the 8 A blocker was 
replaced by a 150A blocker which increased the flux by a factor of nearly 
20 . Again, there was some evidence of a change in angle of ma ximum extinction 
by a degree in angle. Hence, temperature change, if it exists, must be attri- 
buted to a minor IR leakage of the trim filter. In any case, the calcite 
changed by less than . 005 °C or .00l8lA which is certainly not of any conse- 
quence. 

TRANSIENT EFFECTS AND THERMAL GRADIENTS 

By tracking the transmission of the laser beam as a function of time, 
it was possible to measure the temperature of the calcite as a function of 
time. To measure the temperature response on heating the oven was set at 
a control temperature using a temperature controller. The waveplate was 
rotated to obtain a transmission minima, and then a new temperature controller 
was attached to the oven. (Separate temperature controllers were used because 
each individual controller is tuned to control optimally at a specific tem- 
perature. The controller can be fine tuned over a one degree range. ) During 
the heating period, the output of the laser was continuously recorded on a 
strip chart recorder. As the calcite heated, the laser transmission was 
modulated, a new minima occurring when the total phase change of the calcite 
changed tt radians (one fringe). 

When marks on the recorder chart were made at specific temperatures 
indicated by the quartz thermometer during the warm-up, a comparison of 
the thermal response of the filter oven and the calcite crystal was obtained. 
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A plot of oven temperature and calcite temperature versus time during a 
transition from 30° to 36% is shown in Figure 14. As can he seen from 
the curve, there is almost a degree centigrade difference between the aluminum 
block and the center of the calcixe during the heat-up phase. (The laser 
beam probed the calcite at its center. ) For the run shown, the calcixe was 
a press fit into the aluminum housing. 

If the observed temperature difference was occurring across the calcite, 
about 2 fringes would have been observed along the aperture radius. Under 
illumination by a full aperture laser beam, the radial phase change would 
manifest itself by two concentric bright rings which continuously expanded 
out from the center of the crystal as the crystal heated. When illuminated 
by a full aperture collimated beam, no ring structure was observed during 
heat up. To check more quantitatively, the 2mm laser beam monitor was moved 
to pass through the calcite 1. 5 mm from the edge of the crystal and then 
thermal cycle experiment was repeated. The response of the calcite at the edge 
was identical to the response at the center of the crystal. Thus, no thermal 
gradient exists across the calcite even under the fastest heat input condition 
of the controller, .4°C/min. 

To further check the location of the thermal gradient, the element was 
removed and thermal contact silicon grease was installed in the space between 
the crystal, its containing tube, and the oven block. The thermal cycle 
experiment was repeated. The grease reduced the gradient by a factor of two. 
The rest of the gradient most probably is across the alumi num tape that is 
wrapped around the calcite. 

From the tests on the sample element it can be inferred that a- shock 
mounted and space qualified^ calcite element would be isothermal, but that a 
temperature gradient would exist across elastomer mounting substance. The 
size - of the gradient would be dependent upon the thermal input rate. In any 
case, such gradients can be measured and calibrated. 
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TEMPERATURE SENSITIVITY OF CALCITS 


INTRODUCTION 

The sensitivity of the operating wavelength of a birefingent filter to temperature 
is a critical parameter in the design of an active filter controller which corrects 
for temperature changes by tuning elements by rotation of polarizers or wave- 
plates. A knowledge of the temperature sensitivity is also important for the 
design of a constant temperature enclosure. Somewhat surprisingly, the 
temperature sensitivity of calcite has not been measured with high accuracy 
or as a function of wavelength. The only published value being that of Lyot 
at the Sodium D line. 


The thermal effect can be measured by too techniques. The first 
directly measures the change in wavelength of a transmission peak of a Lyot 
element due to a change in temperature. The second measures transmission 
modulation as a function of temperature at a fixed wavelength. The wavelength 
sensitivity is obtained from the transmission modulation through the relation 


where 


_ , 2 2rr 5® P \ 

&T * o bT L / n \ou, 

u a\ ; 


_ o ud 

® T 


(1) 

( 2 ) 


and where p is the index difference, d is the thickness of the crystal, and 
\ the wavelength. The change in phase, <o , with temperature comes from both 
change in index and crystal length. That is. 


6co _ / 1 t 1 ftd (3) 

5T ? % 5T d E,T J 

Comparison of equations (l) and (2) shows that the temperature sensitivity is 
independent of crystal length. 
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Measurements of both the thermal coefficient of expansion and of index are 

given in the International Critical Tables; however, they are not sufficiently 

accurate to determine the wavelength sensitivity. The index values have been 

fitted by a polynomial by Beckers which is accurate to a part in 10 - the 

accuracy of the index values. The temperature derivative term in the polynomial is on 

-5 

the order of 1 x 10 per degree , so that it is not surprising that the 
temperature sensitivity predicted by the critical values are not accurate to 
better than 10-20 percent.- 

To measure the shift in wavelength directly requires a very well calibrated 
high resolution spectrograph. Therefore, these experiments were carried out 
at Kitt Peak. To measure the modulation with temperature requires only stable 
wavelength sources and an accurate photometer. The modulation experiments 
were performed at Lockheed. 

MSASUBEMENT TECHNIQUES 

In order to measure values of b\/b t, a Halle filter oven was modified to 
contain a single calcite crystal. The crystal was a normal wide-field Lyot 
element with a separation between transmission maxima of approximately liL. 

The oven temperature was monitored by the quartz thermometer. The temperature 
was controlled by a set of four Lockheed temperature controllers with nominal 
control temperatures of 30, 36, 40, and 48° C. 

The filter was mounted vertically over the entrance slit of the main specto- 
graph of the McMath Solar Telescope at Kitt Peak National Observatory (which 
is a 13.7 meter Czerny- Turner with a spectral resolution 300,000). The filter 
was fed with an f/60 solar beam. Measurements were made by using the 
spectrograph in a spectrometer mode scanning over approximately 10A centered 
at k-702, 5150, 5850 and 6439A. A total of 12 scans were co-added for each 
measurement. At each of the wavelengths, a solar spectrum was taken before 
and after the measurements at the four temperatures . 

The measurements at Lockheed were made using zhe system used for the thermal 
controller experiments, where the laser could be replaced by a line source 
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with a collimator and appropriate isolation filter (Figure 13). The line source 
measurements were carried out at X7032, 5461, 4800, 4678, and 4358. 

Laser measurements were made at x^328. These measuremenxs used 
the 0.08 a calcite element - the element used in the thermal controller 
experiments . 


KITT PEAK MTIOilAL OESERVATOKI DATA REDUCTION 

All of the raw Kitt Peak data were manipulated using a spectral 
data handling program called REDUCER. The spectral daxa were corrected for the 
grating equation, fourier transformed, filtered, and retransformed. The 
transform operations served to optimally smooth the data. The smoothed data 
were “hen analyzed to find all of the locations of intensity maxima and minima. 
The wavelength differences between corresponding minima or maxima on records 
made at different temperatures yielded 6X/bT- 

The measured values for the temperature sensitivity did not show any temperature 
dependence. The results of the KPNO measurements are shown in Table tv. The 
least squares straight line fit through the data points is 


&X. 

6T 


7.35 x 10' 5 X - .10273 


( 4 ) 


where is in A/° C and X is in angstroms. The mean temperature is 39° C. 

LOCKHEED DATA REDUCTION 

The experiment at Lockheed directly measured the number of intensity cycles between 
two temperatures. Because the number of cycles could be obtained with high 
accuracy it was noted in the first measurements that the fringe count/ 5 C was 
temperature dependent. Since even Lyot had stated that the temperature sensitivity 


B-17 



Table I? 


Temperature Sensixiviry of Calcite Prom KHTO Eata 



4702 

.2433 

5150 

.2771 

5866 

.3306 

6439 

.3704 
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of calcite was independent of temperature from l3°C to 40°C, numerous 
measurements were made at each temperature in order to obtain good statistics. 

Shown in Table V are the mean temperature and mean sensitivity as well as the 
least square fit straight line parameters for each of the measured wavelengths. 

Also shown in Table V are the total number of measurements at each wavelength. 

Using all the straight line fits of the fringe count per degree a~ ~he various tem- 
peratures, a mean -emperature dependent straight line can be obtained for the fringe 
counu. 

F/’C = (1 + 5.88 x 10 ' 8 XT) (5) 

The parameter is excellently fit by 

a - 1.4-94 x 10 '/ \ - .16443 

A. 

Thus the number of fringes/ 0 C is 

F/°C = (-.16^43 + 1.4§4 x 10 V X) (1 + 5-88 x 10 -8 XT) (6) 

where X is in angstroms and T in °C. The values given in equation (6) agree with 
the experimental values to within 0.1% from X 4358 to 17032. The length of the 
calcite element used in the experiment was 13.010 cm, so that the number of 
fringe s/cm/°C is 

F/ cm/° C = (-.01264 + 1.1483 x 10 8 / 1 ) 

x (1 + 5.88 x 10“ 8 T l) (7) 


The fringe sensitivity times the free spectral range (peak to peak separation) 
is just the wavelength sensitivity to temperature. The free spectral range is 


FSR 


1 


2 


Md (1 - 


bM 

bT 



( 8 ) 
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Table V 


Fringe Count Parameters and Temperature Sensitivity 
of Caleite From Lockheed Experiment 



F/°C 

FSR (39°) 

r T < 3sf > 

*0 

II 

cy-^x 10 

N 

4358 

3.31 

.06575 

.2176 



8 

4678 

3.0646 

.07887 

.2405 

3.034 

7.76 

24 

li 800 

2.96 

.08417 

.2492 


- 

8 

5461 

2.602 

.11524 

.2999 

2.568 

8.712 

35 

6328 

2.225 

.1622 

.361 

2.197 

7.114 

60 

7032 

1.994 ' 

.20540 

.4096 

1.958 

9.58 

24 


Note: 1) F/°C = a + a n T 

’ ' o 1 

2) N is the number of measurements 
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The FSR for a centimeter length calculated from the Beckers formula can be 
fitted to better than a tenth of a percent by 

FSH = (-2.722 x 10' 1 - 1.5065 x 10"° X 

+ 5-9725 x 10" 8 x 2 ) (9) 

The product of equations (7) ana (9) represents an analytic expression for the 
temperature sensitivity with an accuracy of a tenth of a percent. 


Thus* ohe final analytic approximation for the temperature sensixivixy is 


= (-3-126 x 10 2 / x + 1.711 x 10~ 3 + 6.86 x 10“ 5 X - 7 .55 x 10“ 10 X 2 ) 

x (1. + 5-88 x 10 -8 XT) (10) 

The temperature sensitivity and the number of fringes/cm/ 3 C versus wavelength 
at 40° are shown in Fig.15. The FSR for a cm of calcite versus wavelength is 
shown in Fig.l6„ 

Although it is not readily apparent from casual comparison of equations (4) and 
(lo)j the KHTO and Lockheed results do agree at 39°C for the region 47OO xo 
65OOA. Both above and below the KENO range (4700-6500) the temperature sensitivity 
departs as can be seen on Figure 15 . 
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AUTOMATIC WAVELENGTH STABILIZATION 


INTRODUCTION 

The traditional source of wavelength drift in birefringent filters has been 
inadequate thermal control of the calcite. The temperature controller 
described above is adequate for control to better than . 005°C for periods 
of weeks to months, which is more than adequate for ground based uses. 

However, for a universal filter with individual tuning controls on each 
element, automatic temperature compensation by tuning is possible. 

Automatic tuning has the advantage of eliminating the temperature controller 
electronics. Absolute wavelength calibration is also possible by reference 
to a laser. 

LOCKHEED UNIVERSAL TECHNIQUES 

The Lockheed Alternate Partial Polarizer Universal Filter (LAFPU) is the 
first filter system designed to operate without a temperature controller. 

Instead, temperature sensors are mounted adjacent to the calcite Crystals 
and the elements are tuned to maintain wavelength based on the signals from 
the sensors and temperature sensitivity of calcite. (Some care has been taken 
in the design to make the filter respond slowly and uniformly to temperature 
perturbations . ) Although additional software complexity is introduced to 
compensate for the lack of a stand-alone temperature controller, the soft- 
ware required is a relatively small addition to wavelength referencing and 
automatic filter assembly which are features which are useful if not necessary 
for a Universal filter. 

For the purpose of auto wavelength referencing, a Helium-Neon laser beam is shown 
■through the filter. The filter control computer can read the values of the laser inten- 
sity transmitted by the filter. This data is used by a control program which successively 
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rotates each element and estimates the position of maximum transmission. Each 
tuning element is set to its individual position of maximum transmission. 

Several iterations of the program are sufficient to optimally assemble the 
filter on \6328. Periodic use of the laser reference and auto- assembly 
guarantees that the filter is operating properly and is on wavelength. 

In order to operate the filter at any other wavelength than X6323 either 
the laser is referenced or the current estimated setting of X 6328 is used 
to set the locations of the waveplates. Since the temperature sensitivity is 
known to a part in 1000, the setting capability is better than a degree in angle 
for the widest element, narrower elements have corresponding higher resolution. 
One degree accuracy corresponds to a wavelength setting accuracy of the 
element FSR/180. 

Operationally, to go to a standard wavelength in the spectrum a table of 
angle offsets from X 6328 is referenced. These values are corrected by the 
temperature state of the filter. That is 


9 k - 9 k + 


180 


FSiV 


l 6328 


bx. 

6 T, 


■6328 


Vk (T . T } 

fseJ ( o) 


( 1 ) 


where $ k (T Q ) is i -ne angular offset from X £>328 of element i- operating at wavelength 
k i a '° uhe ref,erence Temperature T Q , 8* is the angular offset to which element i 
is driven to operate at wavelength k, and FSrJ is the peak to peak separation of 
element i at wavelength k. Any nearby wavelength can be obtained by rotating the 
components proportional to the length of the elements tuned. To maintain wavelength 
the angular positions are updated when a change in temperature is sensed. The 
angular change is 


A 9 - = 


l80 


b T, 


AT 


•k 


( 2 ) 


At present analytic approximations for FSR^_ and exist. Analytic expressions 

T k 

for the now tabulated values of will be developed. They can be easily obtained 

by feeding the filter with a solar beam on the spectrograph bench and using a 
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version, of the wavelength reference auto assembly program that obtains 
intensity data from the spectrometer output. The only requirement is 
that the spectrograph is set in the continuum. 


MEASUREMENTS 

Using the laser ref erence- system built into the Lockheed Universal filter 
each element was successively spun and the position of ma ximum transmission 
measured. After each spin operation, the particular filter tuning element; 
was set to the position of maximum transmission. This operation was carried 
out continously for 6 hours. All spin data was stored. This included the 
set position of all elements, the transmission at the set point, and the 
transmission as a function of angle for each element. After each spin sequence, 
the filter was set to the position of 16328 predicted by the temperature 
compensation program. A comparison of the predicted position and the spin 
set positions yields the relative error of the temperature compensation 
routine. 

During the testing, all stepping motors were always on, even~if not . 
stepping, and the filter was enclosed by a box. In a six hour period the 
interior temperature of the calcite increased from 20°C to 43°C. During the 
entire period the difference between the setting of the filter based on the thermal 
compensation program and the measured setting on the laser line was dominated 
by the least count errors of the thermal sensors. Ignoring least count effects 
the filter stayed on wavelength to within .OIA. for the entire period. 

The test described above is much more severe than ever would be encountered 
in normal operation. If the filter is not enclosed in a box, the interior 
temperature does not rise significantly. Normal operation would not operate 
the motors on much more than a five percent duty cycle. 
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RADIATION EFFECTS 


The radiation problem for birefringent filters comes from two sources. 

One is the incident solar optical radiation and the other is the penetrating 
particle and X-ray flux. 

Because of the sensitivity of the required blocking filters to solar 
radiation below 4-000A and the leakage of polarizers above 7500A, the entrance 
windows must limit the solar optical flux to the 4-000 -7500A band. The filters 
used by the Skylab H-alpha telescope accomplish this and should be used in 
front of the birefringent filter. The solar flux behind the Skylab type filter 
is just the same as normally used on ground-based telescopes. 

Penetrating radiation will have to traverse at least 3 cm of aluminum 
when incident from the side of the filter housing. Radiation along the axis 
of the filter will have all the optics and the secondary mirror structure 
in the path which will be greater than the equivalent of 3 cm of aluminum. 
Measurements of the radiation dose indicates that in a 330 nautical mile polar 
orbit the dosage is well described by 

Dose = 27.4-10 rads/day 

where t is in cm of_ aluminum. Radiation dose in a shuttle orbit would be less 
than this. This relation was obtained by E. Gaines (EMSC) based on data from 
Trapped Particle Environment of SMM , GSFC X- 6 OI- 76 - 97 . For a thickness of 3 cm 
of Al, the dose is 

Dose (3 cm) = 3 x 10 ^ rad/day. 

Radiation damage occurs in polymer, crystals, and glasses at a radiation dose of 

5 T 

ICr to 10 rads ( Space Materials Handbook, Addison-Wesley, 1965 ). 

From the above numbers the likelihood of damage to any of the materials 
inside the housing is very small. 
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SUMMARY - AND CONCLUSIONS 


The "basic direction of this study has been guided by the desire of the 
Quick Reaction and Special Purpose Facility Definition Team for feasibility 
daua on a birefringenx iilter system for an E£x/magnetograph telescope. The urimary 
requirements were for a l/6A or narrower filter that was tunable from Ha to a 
suitable magnetically sensitive line or lines in the visible. The two principal 
unresolved problem areas for a space qualified filter were proper shock mount- 
ing for the calcite filter elements, and wavelength stability. This study 
has developed a prototype tunable filter module which has survived the shuttle 
shake specification. 

Experience has shown that the primary cause for wavelength instability is 
temperature induced wavelength shifts. Tests of the Lockheed Temperature con- 
troller demonstrate sufficient temperature suability for a facility magnetograph. 
Tests of the Lockheed Universal filter using thermal sensors on The calcite 
elements and a software program to control the tuning components have shown 
that it is possible to temperature compensate without a hardware temperature 
controller. An advantage of the software system is the ability to incorporate 
a wavelength reference He-Ne laser to periodically verify the state of the 
filter. The complete software system for wavelength tuning, thermal control, 
and wavelength reference can easily be contained by a small sophisticated desk 
calculator such as an HP 9825. Further studies of the Lockheed Universal filter 
have shown that with a knowledge of the wavelength dependence of temperature 
it is possible to actively tune a filter to any desired wavelength with suffi- 
cient accuracy (> .OlA) without a hardwired temperature control. 

A secondary concern of this study has been the possibility of leakage of 
index marching oil from the filter. An oil sealed container with a novel rotary 
feed-through has been constructed and environmentally tested and demonstrated 
not to leak. 

Oil filled versus non-oil filled filter designs have been evaluated. 

Optical cementing techniques have been developed which allow the majority of 
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the filter components to be mounted in groups, thus minimizing the number 

of air-glass interfaces. Calculation of cemented -A-R coated filter designs 
\ 

have shown that it is practical to eliminate the oil in a birefringent filter. 

A production grease-filled Halle Hex filter has been disassembled and reassembled 
with cemented components and A-R coated surfaces, but without oil. The A-R- 
cemented filter performs satisfactorily. Recently announced A-R coatings 
should allow operation from . 35 to 2 microns without oil. 

The modular approach to filter design and construction illustrated by 
the prototype filter module and the Lockheed Universal filter allows the 
construction of filters of any full width at half -maximum (FWHM) and free 
spectral range (ITSR). Limitations on the full width at half -maximum are 
imposed by length of calcite. However, the existing Lockheed Universal has 
a FWHM a factor of three narrower than that desired by the FDT. The FSR is 
limited only by the number of modules used. The state of the art in blocking 
filters requires six to eight modules for a raagnetograph system. 

As a result of the work performed in this study, in the opinion of the 
authors, a tunable birefringent filter suitable for flight is now within the 
state of the art. However, A-R coated filters are limited to the range of 
'the A-R coating. 
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APPENDIX I 


The purpose of this appendix is to record my (A.M. Title) judgement 
on material and techniques for producing birefringenx elements. 

1. Calcite 

Ax present there seems to he at least four sources of calcite. The 
calcite comes from mines in Mexico, South Africa, India, and perhaps the USA. 

We have conxacxed a half doaen agents and have obtained prices in the range 
of $750 to $2000 per kilogram. The calcite in the Lockheed Universal, as well as 
all xhe Spectra Optics filters, has come from Mexican mines represented by Gary 
Anderson. Karl Lambrect Company sells both Mexican and South African calcite. 

We have obtained a sample of the South African material and it is of signi- 
ficantly higher whiteness of any previous Mexican calcite. Lambrect claims 
to have a significanx stock of such material and guarantees its quality 
at $2000 1 kg. 

Harry Ramsey has had extensive conversaxions with representatives or 
owners of calcite mines and sample calcite has been examined. If high quality 
calcite exists, there is no reason to look further for alternative materials. 

It would appear prudent therefore to have someone knowledgeable about calcite 
and calcite working actually visit and examine to known mined stock piles 
to assess the claims of xhe mine owners and their representatives. Most 
critical is the sxock of mined Souxh African calcite, since xhis material, 
because of political problems, could rapidly become unavailable. In view 
of the high quality of the samples in hand, it is my opinion that sufficient 
calcite can be obtained. 

2. Metal Acid Pthalates 

The most attractive alternate crystal appears to be potassium acid 
pthalate (KAP) or one of its alkali metal relatives. These crystals have been 
grown in fairly large sizes and have an effective birefringence nearly thax 
of calcite. We have some of the crystals produced at Hew Mexico State 
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and round them to be of 1 non-uniform birefringence. These crystals are water 
soluable, fragile, and extremely sensitive to therm al shock. 

At present, I feel that investigating KAP and ins relatives is a lower 
priority nhan the plastic materials discussed below. 

3. Plastic Materials 

Polyvinyl alcohol ( PVA) is a major commercial plastic. It is the plastic 
film base material for Polaroid Corporation polarizer and waveplate sheets. 

PVA becomes birefringent when snrenched. Reasonable snretches yield birefrin- 
gencies up to three times quartz. The plastic can be stretched very uniformly. 
The sheen material can be la m inated to form birefringent elements. We have 
made a 80A tunable filter using FVA sheet. The sheet and the laminate are scanner 
free and are obtainable in at least 15 cm diameter. 

FVA sheen should be considered for the purpose of making tunable pre- 
filters for the final filter. Sole type filters become much more attractive 
using laminate sheets. I would recommend looking inno the Universal blocking 
filter using plastic sheet. 

Michelson Elements as Analog Lyot Elements 

Evans in 1949 suggested that Michelson interferometers could be used 
to build Lyot filters. The response of a fixed Michelson is identical to 

that of a nyot element. A fixed wide field Michelson offers two significant 
advantages for elements narrower than 1 /4 A. First, since a Michelson is 
roughly cubical regardless of 'its effective FW5M for elements narrower than 
l/b A, there is increasing greater length advantages. Second, a wide field 
Michelson has a field of view which is the square of the field of an 
identical bandpass wide field Lyot element. 

We have built a prototype l/20 A polarizing wide field Michelson. The 
device is quite stable and easily tuned by rotation of a waveplate or polarizer, 
a 1/50 A filter using 2 or 3 wide field Michelsons is a reasonably straight- 
forward development. 

m my opinion j a l/50 A tunable filter could almost eliminate the need 
for high dispersion spectrographs. Xn would provide an enormously powerful 

tool for line profile analysis and magnetic field strength. 
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FIGURE CAPTIOUS 


Figure I 
Figure 2 
Figure 3 

Figure k 

Figure 5 

Figure 6 
Figure 7 
Figure 8 

Figure 9 
Figure 10 

Figure 11 


Figure 12 

Figure 13 
Figure 14 

Figure 15 
Figure 16 


Pallet random vibration levels. 

Filler random vibration levels. 

Fringe patterns visible in preliminary test unit due to separa- 
tion of waveplates. 

Mechanical drawing of Mark I and II Preliminary Test Unit (dimen- 
sions in cm) . 

Photograph of Mark II Preliminary Test Unit and mounting fixture 
for shake table. 

Mechanical drawing of prototype filter module (dimensions in cm). 
Photograph of prototype filter module prior to assembly. 

Exploded photograph of prototype filter module and sealed 
enclosure. 

Photograph of assembled sealed enclosure. 

Transmission versus number of modules for several values of 
reflection loss per surface. 

Reflection coefficient (full scale one percent) versus wavelength 
for a simple four layer A-R coating. Properties calculated using 
a published (Thetford, Optica Acta 11 , 159 , 1969 ) design and 
program thin film. 

Drawing of Halle oven showing location of the quartz probe and 
controller sensors. 

Optical schematic for the thermal controller test system. 

Temperature response of the quartz probe (oven) and calcite interior 
versus time. 

Temperature sensitivity and fringes/cm/°C versus wavelength at 40°C, 
Free spectral range versus wavelength for a one centimeter calcite 
element . 
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